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1.  INTRODUCTION 


BRLCB  is  a  comprehensive  closed-chamber  data  analysis  program  designed  to  perform  all  analysis 
associated  with  the  high-pressure  combustion  research  facility  located  at  the  Weapons  Technology 
Directorate,  U.S.  Army  Research  Laboratory  (ARL).  The  program  offers  five  different  analysis  modes. 
First,  the  traditional  bum  rate  analysis  which  determines  the  apparent  linear  burning  or  regression  rate  for 
a  propellant  given  an  experimental  closed-chamber  pressure  .iistory  together  with  the  propellant 
thermochemistry  and  geometry.  Second,  a  synthetic  prcssure-timf;  generation  mode  which  calculates  the 
expected  closed-chamber  pressure- time  profile  given  propellant  thermochemistry,  linear  bum  rate,  and 
geometry.  The  third  analysis  mode  is  determination  of  propellant  surface  area.  This  analysis  dctcmiines 
the  reacting  (burning)  surface  area  of  the  proprcllant  required  to  support  the  experimental  pressure-time 
profile  given  propellant  thermochemi.stry,  geometry,  and  linear  burning  rate.  Analysis  mode  4  is  referred 
to  as  an  interrupted  chamber  analysis  which  is  identical  to  the  bum  rate  analysis  except  that  the  pressure- 
time  profile  is  for  incomplete  (interrupted)  prop)cllant  burning  which  results  when  a  "blow-out  disc"  in  the 
closed  vessel  ruptures  at  a  predefined  pressure  and  terminates  propellant  combustion.  Tire  fifth  and  final 
mode  of  analysis  is  also  a  bum  rate  analysis,  but  is  for  the  analysis  of  propellant  burning  coupled  with 
the  injection  of  an  electrically  generated  plasma  into  the  closed  chamber.  This  mode  is  specifically 
designed  for  investigation  of  potential  propellants  for  use  in  clcchothemial-chemical  (ETC)  gun 
application,'.  In  addition  to  the  five  analysis  modes,  the  program  includes  extensive  provisions  for  uata 
preparation  (smoothing,  etc.),  graphics  capabilities,  and  pwst-proccssing  for  prcparatioii  of  output. 

Over  the  years,  a  number  of  programs  have  been  dcvclop)cd  to  perform  both  the  bum  rate  and  surface 
area  analysis  (Robbins  and  Horst  1976;  Price  and  Juhas/.  1977;  Obcrlc,  Juhasz,  and  Griffic  1987;  Obcrle 
and  Kookcr  1989).  However,  except  for  the  Obcrle  and  Kookcr  program,  lliesc  programs  were  designed 
to  ,.ccommodatc  oniy  homogenous  proptcliants  (i.c.,  propellants  with  constant  thcnnochemical  properties). 
An  impertant  feature  of  BRLCB  is  die  ability  to  analyze  not  only  homogeneous  propellant  burning  but 
also  layered  and  deterred  pmpeilant  burning. 


The  principal  objective  of  this  repert  is  to  provide  details  for  u.sing  the  various  program  options  and 
analysis  modes.  In  addition,  the  rcpert  provides  details  of  the  theory  underlying  the  different  analysis 
modes,  results  of  test  ca.scs  utilized  in  verifying  and  lc.sting  the  computer  program  implementation  of  the 
model,  and,  finally,  comparison  with  results  obtained  using  other  closed-chamber  data  analysis  programs. 
To  accomplish  this  objective,  this  rcjxnl  is  divided  into  four  mam  sections;  (1)  User’s  Guide  and  Program 
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Description,  (2)  Theoretical  Analysis  of  Gosed-Chamber  Combustion,  (3)  Validation,  and  (4)  Comparison 
With  Other  Gosed-Chamber  Data  Reduction  Programs. 


2.  USER’S  GUIDE  AND  PROGRAM  DESCRIPTION 

2.1  Program  Structure  and  Overview.  BRLCB  is  intended  for  use  on  an  IBM  PC  or  compatible. 
However,  the  program  is  written  in  standard  FORTOAN  77  and  can  be  ported  to  operate  on  any  system 
which  supports  FORTRAN,  the  exceptions  being  the  graphics  routines  utilized  in  the  program.  In 
actuality,  the  BRLCB  program  structure  is  not  that  of  single  large  program  but  of  8  individual  FORTRAN 
programs,  1  DOS  batch  file  program,  and  34  ancillary  information  files  which  are  listed  in  Table  1. 
(Note:  To  run  BRLCB,  all  of  the  files  listed  below  should  be  contained  in  a  separate  directory,  preferably 
on  a  hard  disk.  In  this  report,  it  will  be  assumed  that  all  the  files  arc  in  a  separate  directory  with  the  user 
logged  into  that  directory.) 

Table  1.  Required  FORTRAN,  DOS  Batch,  and  Ancillaiy  Files  for  BRLCB 


FORTRAN  Programs 

MKCHCE.EXE 

MKMASTER.EXE 

MKGAGE.EXE 

MKPl’DATA.EXE 

MfJNF-EXE 

MKSMOOTH.EXE 

MKCAL.EXE 

MXOUT.EXE 

DOS  Batch  File 

BRL.BAT 

Ancillary  Files 

IGNITER.INF 

GAGEFILE 

S5.CFF 

D5.CFF 

S7  CFF 

D7,CFF 

S9.CFF 

D9.CFF 

Sll.CFF 

Dll.CFF 

S13.CFF 

D13.CFF 

S15.CFF 

D1.5.CFF 

S17.CFT 

D17.CFF 

S19.CFF 

D19.CFF 

S2LCFF 

D21.CFF 

S23.CFF 

D23.CFF 

S25.CFF 

D25.CFF 

S27.CFF 

D27.CFF 

S29.CFF 

D29.CFF 

S3 1. CFF 

r3LCFF 

S33.CFP 

D33.CFF 

S35.CFF 

D35.CFF 

A  complete  listing  of  all  the  programs  and  files  arc  found  in  Appendices  D-L.  A  brief  description 
of  each  program  or  file  is  provided  below. 

MKCHCE.EXE  :  Overall  controller  program  for  BRLCB  which  allows  program  options  to  be 
invoked  from  a  menu. 

MKMASTER.EXE  :  Program  used  to  create  master  information  file  to  be  used  in  the  data  analysis. 

MKGAGE.EXE  :  Program  used  to  maintain  pressure  transducer  calibration  infomiation  file 
which  is  used  if  voltage-time  data  must  be  converted  to  prcssure-time  data. 

MKPTDATA.EXE  :  Program  used  to  import  experimental  data  into  BRLCB. 

MK1NF.EXE  :  Program  used  to  update  a  master  information  file  created  in  program 
MKMASTER.EXE  for  a  specific  closed-chamber  experiment. 

MKSMOOTH.EXE  :  Program  used  to  prepare  the  pressure-time  data  for  the  analysis. 

MKCAL.EXE  :  Program  which  performs  the  selected  data  analysis. 

MKOUT.EXE  :  Program  which  provides  both  paper  and  disk  output  of  the  program  rcsulLs. 

BRL.BAT  :  DOS  batch  file  used  with  MKCHCE.EXE  to  control  program  flow. 

IGNITER. INF  ;  Data  file  which  stores  thcrmochcmical  properties  of  various  igniters. 

GAGEFILE  :  Data  file  which  stores  pressure  transdutujr  calibration  information. 

S5.CFF — S35.CFF  :  Data  files  containing  the  numerical  cocfficierts  utilized  in  smoothing  the 
pressure-time  data. 

D5.CFF — D35.CFF  :  Data  files  containing  tlic  numerical  coefficients  utilized  in  differentiating  the 
pressure-time  data. 

The  remainder  of  this  section  will  cover  running  the  program  and  options  which  must  be  run  for  the 
different  analysis  modes.  Subsequent  sections  will  provide  additional  details  on  the  seven  main  programs 
which  comorisc  BRLCB. 

To  begin  BRLCB,  at  the  DOS  prompt  type  BRL.  This  will  invoke  the  batch  file  listed  above,  which 
after  some  initialization  will  call  the  program  MKCHCE.EXE.  The  main  BRLCB  menu  (sec  Table  2)  will 
tlien  be  displayed. 


Table  2.  Main  Menu  for  BRLCB 


1  BRL.CB  **  Version  3.0  •*  January  1992 

Main  Menu 

1. 

Create  Master  Information  File 

2. 

Update  gage  information 

3. 

Prepare  pressure-time  data 

4. 

Prepare  firing  information  file 

5. 

Smooth  pressure-time  data 

6. 

Perform  data  analysis 

7. 

Prepare  output 

8. 

Exit 

Please  enter  your  choice  (1-8). 

The  particular  sequence  of  options  to  select  will  depend  upon  the  analysis  mode  desired.  Figures  1 
and  2  illustrate  the  option  flow  for  the  five  analysis  modes.  The  options  must  be  performed  in  the 
indicated  order  if  the  program  is  to  work  correctly,  and  all  options  indicated  must  be  selected.* 
Appendix  A  lists  the  information  required  by  tlic  user  for  each  option. 


*Option  1.  Crcaw  master  information  file,  docs  not  have  to  be  invoked  if  an  aiipropriate  master  file  for  the  analysis  already 
exists. 


Option  1 

Create  Master  File 


Figure  1.  Oplior.  flow  for  analysis  dodes:  (1)  bum  rate  reduction:  (3)  surface  area  anal 
interruDted  burner  and  (5)  ETTC  bum  rate  reduction. 


igurc  2.  Option  flow  for  analysis  mode  (2)  generate  pressure-tiine  data. 


2.2  Option  1,  Create  Master  Information  File  (T>rogram  MKMASTCR.EXE).  In  BRLCB,  the  master 
information  file  is  the  starting  point  regardless  of  the  analysis  mode.  The  idea  behind  the  master 
information  file  (.MAS  suggested  file  extension)  is  tliat  certain  infomiation  (e.g.,  propellant  geometry)  will 
remain  fixed  for  several  experimental  firings.  Thus,  one  master  file  will  be  created  which  will  then  be 
utilized  to  create  an  updated  file  (referred  to  as  the  .INF  in  Option  4,  Preparing  Firing  Information  file). 
The  .INT  file  will  contain  information  specific  (e.g.,  propellant  mass)  to  a  single  experimental  firing. 

Figure  3  provides  an  outline  of  the  information  required  in  Option  1. 


MAKMASTER.EXE 

- IdentIUcatlon  Information 


—  Propellant  Information 

A.  Identification 

B.  Thermochemistry 

1.  Number  of  Layers 

2.  Constant  Properties  vs.  Varying  Properties 

Flame  Temperature  (I<) 
impetus  (J  /  g) 

Molecular  Weight 
Covolume  (cc  /  g) 

Ratio  of  Specific  Heat  (-) 

Density  (g  /  cc) 

C.  Propellant  Geometry 


Igniter  Information 


Write  Master  File  (.MAS  suggested  file  extension) 


Bxit 


Figure  3.  Overview  of  information  required  in  Option  I. 


Table  3  displays  the  first  screen  presented  when  Option  1  is  run.  As  shown  in  the  figure,  llic  master 
information  file  can  be  created  from  an  existing  master  information  file,  However,  this  capability  is  very 
limited — editing  single  entries  of  an  existing  rna.sicr  information  file  —is  not  available  nor  can  the 
infomiation  contained  in  the  previous  master  infomiation  file  be  previewed.  In  building  from  an  existing 
file,  it  is  assumed  tliat  the  user  is  cognizant  of  the  contents  of  the  previous  master  information  file  and 


Tabic  3.  Option  1.  Screen  1 


BRLCB  **  Vereion  3.0  **  January  1992 

Creating  a  Master  Information  file 

This  program  will  create  a  master  information  file  which 
will  contain  information  that  does  not  change  from  one 
firing  to  the  next. 

Create  .MAS  file  from  existing  file?  [Y/N] 

Default:  Y 


will  use  the  file  to  avoid  rc-cnicring  certain  blocks  of  information  (generally,  thermochemistry  or  geometry 
under  propellant  infomiation,  Figure  3).  It  is  recommended  that  master  information  files  be  created  from 
the  beginning  until  the  user  becomes  very  familiar  with  the  overall  BRLCB  program  structure.  The 
structure  of  master  informatinn  and  .INF  (created  in  Option  4)  files  is  provided  in  Appendix  B.  These 
files  arc  saved  as  ASCII  files  and  could  be  directly  edited  with  any  standard  word  processor,  bypassing 
Options  1  and  4  altogether.  However,  care  must  be  taken  in  editing  tlicsc  files  since  all  subsequent  file 
names  (automatically  generated  by  the  program)  utilized  in  die  analysis  arc  stored  in  these  files.  Failure 
to  correctly  modify  any  of  the  stored  file  names  will  result  in  loss  of  data.  Direct  editing  of  the  master 
information  or  .INF  file  is  not  recommended. 

2.2.1  Identification  Information.  After  selecting  whether  or  not  to  create  the  master  information  file 
from  an  existing  file,  the  identification  information  (Figure  3)  suboption  is  automatically  selected.  It  is 
in  this  suboption  that  ilic  file  name  associated  wiili  tlic  master  information  file  being  created  is  requested. 
Although  any  acceptable  DOS  name  and  extension  can  be  u.scd,  IT  IS  RECOMMENDED  THAT  ALL 
MASTER  INFORMATION  FILES  HAVE  AN  EXTENSION  OF  .MAS.  Once  the  idenUfication 
information  is  completed,  the  user  is  returned  to  the  main  menu,  as  shown  in  Table  4. 

2.2.2  Propellant  Infoimation  -  Identification.  As  shown  in  Figure  3,  the  propellant  information 
suboption  has  three  parts — identification,  thermochemistry,  and  geometry.  Identification  requests 
information  concerning  the  propellant  type,  source,  and  lot.  A  complete  understanding  of  the  information 
requested  and  tlic  choices  m.ide  in  the  ihcnnochcmistry  and  geometry  sections  of  tlic  propellant 
information  suboption  arc  required  for  effective  use  of  BRLCB.  A  detailed  description  for  these  sections 
is  provided  in  the  next  two  sections  of  this  report. 
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Tabic  4.  Main  Menu  for  Option  1 


Creation  of  Master  File:  Main  Menu 

1. 

Identification  information  * 

2. 

Propellant  information 

3. 

Igniter  information 

4. 

Write  ilic  master  file 

5. 

Exit 

An  * 

indicates  tliat  the  information  has  been  provided. 

Please  enter  your  choice  (1-5). 

(1)  Propellant  Information  -  Thermochemistry.  In  tiiis  poilion  of  the  progiani,  thcnnochemical 
properties  associated  with  a  single  grain  arc  supplied  to  the  program.  BRLCB  makes  no  distinction 
between  homogeneous  (single  layer  with  constant  tiicrmochcmical  properties),  layered  (several  layers  witli 
each  layer  having  constant  thermochemical  properties),  and  deterred  (continuously  varying  thcrmochcmical 
properties  in  at  least  one  layer)  propellant  grains.  However,  the  required  computational  time  for 
homogeneous  and  layered  propellants  is  several  orders  of  magnitude  faster  than  for  variable 
thcnnochemical  properties,  as  shown  in  Table  5,  where  a  five-layer  sphere  with  constant  properties  in  each 
layer  is  analyzed  as  a  layered  propellant  and  then  as  a  deterred  propellant.  Thus,  whenever  possible, 
deterred  propellants  should  be  described  as  a  multilayered  propellant.  Describing  a  propellant  as 
consisting  of  15  layers  with  constant  thcmiochemical  pro[>crtics  in  each  layer  will  result  in  a  substantially 
faster  run  time  than  describing  the  propellant  grain  as  a  single  layer  with  continuously  varying 
thcrmochcmical  properties. 

Table  5.  Run  Times  on  CRAY  XMP  for  Five-Layer  Sphere,  Constant 
Thcrmoclicmical  Properties  in  F,ach  Layer 


Method 

CPU  Time 

Layered 

0.206 

Deterred 

31.6 

The  first  information  requested  under  the  ihcnnochcmisiry  .section  is  for  tire  number  of  layers  in  the 
grain.  Tlic  grain  can  have  up  to  15  layers.  If  the  number  of  layers  is  greater  than  one,  tlie  user  is 
prompted  to  enter  the  starting  depth  of  each  layer,  as  shown  in  Table  6.  Note  from  Table  6  that  grain 
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slivering  should  only  occur  in  the  inner  most  layer  of  the  grain.  However,  BRLCB  will  still  correctly 
perform  the  desired  analysis  if  grain  slivering  occurs  in  a  layer  besides  the  innermost  layer.  The  authors 
feel  that  determination  of  layer  depths  and  changing  thermochemical  properties  when  slivering  occurs  (in 
actual  propellant  grains)  will  be  very  inaccurate. 


Table  6.  Prompt  for  Starling  Depth  of  Each  Layer 


Enter  the  starting  depth  for  each  layer.  The  first  layer 
starts  at  a  depth  of  0  cm  and  will  be  automatically  entered. 
The  last  layer  must  start  at  a  depth  no  deeper  tliat  one-half 
the  length  of  the  smallest  web  (i.e.,  no  slivering  may  occur 
except  in  the  inner  layer  of  the  grain). 


To  illustrate  how  BRLCB  interprets  layer  starting  depth,  consider  Figures  4  and  5.  Figure  4  is  to 
represent  propellant  grains  which  have  no  perforations  or  slots  (in  this  case,  a  spherical  grain  of  radius 
9  cm  consisting  of  three  layers).  Layer  1  is  always  the  outside  layer  and  .starts  at  a  depth  of  0.0  cm. 
Successively  numbered  layers  proceed  into  the  grain  with  the  highest  numbeied  layer  being  the  innermost 
layer.  As  shown  in  Figure  4,  the  starting  depth  for  layer  two  is  4.0  cm  and  for  layer  three,  7.0  cm. 
Figure  5  is  to  represent  propellant  grains  with  perforations.  Specifically,  Figure  5  shows  the  end  view  of 
a  single,  perforated  grain  with  three  layers.  Note  that  layers  1  and  2  form  two  distinct  regions  into  the 
grain.  Layer  1  is  tlic  layer  adjacent  to  the  propellant  grain  boundary  (in  this  case,  the  perforation 
boundary  and  outer  boundary),  while  layer  3  is  the  inneimost  layer.  Layer  1  starts  at  a  depth  of  0.0  cm, 
layer  2  at  1.5  cm,  and  layer  3  at  2.5  cm. 

After  the  number  of  layers  and  starting  depth  of  each  layer  is  entered,  the  user  is  prompted  to  enter 
whether  the  layers  will  all  have  constant  thcrmochcmical  properties  or  whether  at  least  one  layer  will  have 
continuously  varying  properties.  If  all  layers  have  constant  ihermochemical  properties,  the  values  for  the 
properties  listed  in  Table  7  will  be  requested  for  each  layer. 
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starting  Depth  s:  0,0  cm. 


Layer  2 

Starting  Depth 


Layer  3 

Starting  Depth 


1.6  cm. 


2.5  cm. 


Perforation 


Table  7.  Required  Thcrmochemical  Properties 


Property 

Units 

Flame  temperature 

K 

Impetus 

J/g 

Molecular  weight 

- 

Covolume 

cmVg 

Ratio  of  specific  heats 

- 

Density 

g/cm^ 

A  value  for  molecular  weight  consistent  with  the  flame  temperature  and  impetus  will  be  suggested  by 
the  program.  This  value  should  be  accepted,  it  is  the  flame  temperature  and  molecular  weight  which  is 
used  in  the  analysis,  not  the  impetus.  Tlic  molecular  weight  is  computed  from  the  flame  temperature  and 
impetus  by  the  following  equation 

I,,  .  (1) 

Impetus 

If  any  layer  has  continuously  vaiying  properties,  the  user  will  be  prompted  (sec  Table  8)  to  enter  the 
thcrmochemical  properties  shown  in  Table  7  at  both  the  beginning  and  end  of  each  layer.  The  molecular 
weight,  in  this  ease,  is  not  computed  by  the  program:  the  user  must  determine  a  consistent  value  with 
flame  temperature  and  impetus  for  the  molecular  weight.  Equation  (1)  can  be  utilized  for  this  purpose. 


Table  8.  Prompt  for  Entering  Continuously  Varying  Thcmtiochcmical  Properties 


Values  will  be  entered  for  the  beginning  and  end  of  each  layer.  All  input 
in  metric  units.  For  any  layer,  there  should  be  no  more  than  a  20% 
variation  in  any  property. 


(2)  Propellant  Information  -  Geometry:  Currently,  BKLt„i>  supports  13  grain  geometries  ,  as  showm 
in  Table  9. 
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Tabic  9.  Grain  Geometries  Supported  by  BR1,CB 


1 .  Sphere 

2.  Coni 

3.  Rectangular  strip 

4.  1 -Perforated  cylinder 

5.  Slotted  tube 

6.  7-Perforatjon  cylinder 

7.  7-Perforation  hexagonal 

8.  19-Perforation  cylinder 

9.  19-Perforation  hexagonal 

10.  37-Pcrforation  hexagonal 

1 1 .  Cord  with  inhibited  ends 

12.  Sandwich  with  inhibited  sides 

13.  Cigarette 


Implementation  of  the  first  10  grain  geometries  is  identical  to  that  found  in  the  interior  ballistic  code 
IBHVG2  (Anderson  and  Fickie  1987).  Figures  6-18  show  each  of  the  grain  geometries.  (Figures  6-15 
are  from  the  Anderson  and  Fickie  report.)  All  perforations  on  a  grain  are  assumed  to  have  the  same 
diameter.  For  hexagonal  grains,  the  inner  and  middle  webs  must  be  equal.  The  following  notation  is  used 
in  the  diagrams. 


D 

;  Grain  diameter 

WI 

:  Inner  web 

L 

;  Grain  length 

WM 

:  Middle  web 

DP 

;  Perforation  diameter 

WO 

;  Outer  web 

Once  the  propellant  grain  dimensions  are  entered  by  the  user,  the  program  will  determine  if  the 
dimensions  are  consistent.  If  not,  the  program  will  provide  the  user  with  a  consistent  set  of  dimensions. 
The  user  has  the  option  to  accept  the  consistent  set  or  keep  the  dimensions  entered.  However,  for 
computational  purposes,  the  dimensions  must  be  consistent  and  the  program  will  utilize  the  consistent  set 
of  dimensions  suggested  to  the  user  in  the  computation  even  if  the  user  elected  to  keep  the  inconsistent 
dimensions. 
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Figure  14.  19-Pcrforation  Hexagonal.  Figure  15.  37-Pci1oraiii)n  Hexagonal. 


Figure  17,  Sandwicli  With  Inhibited  Sides. 


Figure  16.  Cord  With  Inhibited  Ends. 


Burning  Only 
on  End  Surface 


Figure  18.  Ciiiarctte. 

2.2.3  Igniter  Information.  Three  options  for  igniter  information  arc  provided  by  BRLCB,  as  .shown 
in  Table  10. 


Table  10.  BRLCB  Igniter  Options 


Creation  of  Ma.stcr  File:  Igniter  Information 

Information  concerning  the  type  of  igniter  being 
used  is  needed.  The  available  options  arc: 

1.  Black  powder 

2.  Enter  new  information 

3.  Use  igniter  library 

4.  Exit 

Please  enter  your  choice. 


Option  1,  Black  Powder,  sets  the  igniter  to  black  powder  with  ihcrmochcmical  properties  shown  in 
Table  11.  Selecting  "Enter  New  Information"  will  require  the  user  to  enter  the  same  ihcnnochcmical 
values  as  for  the  propellant  shown  in  Table  7.  ''U.sing  Igniter  Library"  acccs.scs  the  igniter  information 
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file  where  information  on  various  igniters  has  been  saved.  (Note:  ITie  information  in  the  igniter 
information  file  is  built  up  by  the  user.) 


Table  11.  Thermochcmical  Properties  of  Black  Powder  Utilized  in  BRLCB 


Theo.  impetus  (J/g): 

290.00000 

Flame  temperature  (K): 

2188.00000 

Density  (g/cm^); 

1.75000 

Average  molecular  weight: 

66.37000 

Covolumc  (cm^/g): 

0.78500 

Gamma: 

1.21840 

After  the  igniter  information  is  entered,  the  program  returns  to  the  main  menu,  Table  4,  Option  4, 
Write  Master  File  should  then  be  selected.  [Note:  To  preview  tlie  thermochcmical  properties  entered  prior 
to  writing  the  Master  File,  enter  a  10  for  your  selection.) 

2.3  Option  2,  Update  Gaec  Information  (Prottram  MKGAGE.EXE).  Generally,  closed-chamber 
pressure  transducers  record  "counts"  or  voltage  which  must  be  converted  to  units  of  pressure,  which  for 
this  program  must  be  mcg;tpascals  (MPa).  BRLCB  permits  data  to  be  imported  to  the  program  either  in 
units  of  pressure,  voltage  units  (or  any  units  which  can  be  converted  to  pressure),  or  in  a  format  specific 
to  BRL  (details  provided  in  section  on  Option  3,  Prepare  Pressure-Time  Data).  However,  since  each 
research  facility  will  have  their  own  standard  operating  procedure  for  data  acquisition,  only  the  simplest 
(except  for  the  BRL  procedure)  voltage -to-pressure  conversion  is  provided  in  BRLCB.  This  conversion 
is  based  upon  a  second-order  relation  between  voltage  and  pressure.  Tire  conversion  is  given  by 


Pressure  (MPa)  =  A  +  Bx  +  Cx'^ 


(2) 


where  X  is  whatever  quantity  the  user  is  employing  to  import  the  data.  The  coefficients  A,  B,  and  C  will 
depend  upon  the  units  of  the  quantity  x.  However,  the  pressure  must  have  units  of  megapascals. 

The  purpose  of  Option  2,  Figure  1,  is  to  provide  a  utility  by  which  the  user  can  maintain  a  data  file 
of  tlic  conversion  coefficients  A,  B,  and  C,  which  at  the  Weapons  Technology  Directorate,  ARL,  are 
related  to  tlic  pressure  gages.  ITiis  data  file  can  titen  be  accessed  in  Option  3  (Prepare  Pressure-Time 


18 


Data)  when  the  data  are  imported  to  BRLCB.  Using  this  data  file  will  minimize  the  number  of  times  the 
coefficients  must  be  entered  into  the  program  if  the  same  coefficients  will  be  used  for  several  data  sets. 
Maintenance  options  available  for  manipulating  the  gage  file  arc  presented  in  Table  12.  These  same 
options  arc  also  accessible  in  Option  3  (Prepare  Pressure-Time  Data). 


Table  12.  Gage  Maintenance  Options 


G.igc  Maintenance  Program 

Sub  Menu 

1. 

View  gage  information 

2. 

Add  a  gage  to  data  base 

3. 

Delete  a  gage  from  data  base 

4. 

Locate  specific  gage  in  data  base 

5. 

Exit 

Please  enter  your  choice  (1-5). 

2.4  Option  3.  Prepare  Prc.ssurc-Time  Data  (Protiram  MKPTDATA.EXE).  As  indicated  in  the  previous 
section,  BRLCB  provides  five  methods  by  which  data  can  be  imported  to  tlic  program.  These  five 
methods  arc  shown  in  Table  13. 


Table  13,  Methods  by  Which  Data  Can  Be  Imported  to  BRLCB 


BRLCB  supports  five  options  for  preparing  the  pressure-time  data  for 
the  compulation. 

1.  ASCII  file  of  time  and  pre.ssurc  (two  columns,  time/pressurc) 

2.  ASCII  file  of  pressure  (one  column,  pressure) 

3.  ASCII  file  of  lime  and  voltage  (two  columns,  Ume/voltagc) 

4.  ASCII  file  of  voltage  (one  column,  voltage) 

5.  A  voltage-lime  file  from  VuPoint/BRL  procedure 

6.  Exit  option 

How  will  the  pressure-time  data  be  entered? 

Please  enter  your  choice  (1-5  or  6  to  EXIT), 


After  tile  method  by  which  the  data  will  be  imported  to  the  program  is  selected,  the  file  name  for  the 
data  is  requested.  This  will  be  followed  by  a  prompt  for  the  file  name  under  which  the  converted  data 
will  be  saved.  The  prompt  is  shown  in  Table  14,  As  stated  in  the  prompt  (Table  14),  the  file  is  to  have 


no  extension  (for  example,  aitestl,  an  extension  of  .pvt  will  be  automatically  appended  by  the  program 
to  create  the  file  aitcstl.pvt).  The  format  of  the  .pvt  file  is  shown  in  Table  15.  This  file  name,  in  the 
example  a:testl,  will  be  the  file  name  used  with  different  extensions  (sec  Section  2.5,  Option  4)  for  all 
the  files  created  during  the  analysis.  IN  ALL  SUBSEQUENT  PROGRAM  OPTIONS  WHEN  A  FILE  NAME 
IS  REQUESTED,  THIS  FILE  NAME  (WITHOUT  EXTENSION)  SHOULD  BE  USED.  THE  ONLY 
EXCEPTION  BEING  THE  MASTER  FILE  NAME  REQUEST  IN  OPTION  4. 


Table  14.  File  Name  Prompt 


Table  15.  Format  of  .PVT  File 


Line  1: 

Number  of  data  points 

Line  2: 

Maximum  pressure  (MPa) 

Line  3: 

Minimum  pressure  (MPa) 

Line  4: 

Ending  time  (s) 

Line  5; 

Starting  time  (0.0  s) 

Line  6: 

5  -t-  number  of  data  points  N: 

Col  1 ;  time(s).  Col  2:  pressure  (MPa) 

Line  N-t-1; 

Gage  ID  ("none"  if  no  gage  used) 

Line  N+2: 

Voltage  input  to  gage  (0.0  if  no  gage  used) 

Line  N+3: 

Constant  coefficient  (0.0  if  no  gage  used) 

Line  N+4: 

Linear  coefficient  (0.0  if  no  gage  used) 

Line  N-r5; 

Second-order  coefficient  (0.0  if  no  gage  used) 

Once  the  file  names  have  '  ecn  entered,  a  message  describing  the  format  of  the  imported  data  file  is 
provided.  (Note:  For  method  5,  a  voltage-time  trie  from  VuPoini/BRL.  Procedure,  see  Section  2.4.1).  For 
all  other  methods  of  importing  data,  the  total  number  of  data  points  cannot  exceed  2,048  points.  Next, 
Utc  time  step  utilized  in  the  data  recording  is  requested.  This  time  step  is  to  be  given  in  milliseconds  (rns) 
and  the  data  must  have  been  recorded  with  this  constant  time  .step.  The  burning  rate  and  .surface  area 
analysis  are  predicated  upon  a  constant  time  step  for  the  data  recording.  This  program  is  not  de.signed 
to  analyze  data  with  unequal  time  steps. 
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If  the  imported  data  arc  already  in  pressure  units,  the  output  file  will  be  written  and  the  option  will 
be  completed.  For  imported  data  in  voltage  (or  any  other  units),  conversion  coefficients  for  a  second- 
degree  relation  will  be  requested.  These  coefficients  can  be  entered  directly  into  the  program  or  the  gage 
file  (see  section  2.3)  can  be  accessed  to  obtain  tlic  coefficients.  Provisions  for  saving  newly  entered 
coefficients  into  the  gage  file  arc  provided.  Once  the  coefficients  are  entered,  the  output  file  will  be 
written  and  the  option  will  be  completed. 

2.4.1  VuPoint/BRL  Procedure.  This  method  is  based  upon  tlie  data  recording  procedure  utili/ed  at 
the  Weapons  Technology  Directorate,  ARL.  Data  are  recorded  on  Nicolet  oscilloscopes  and  then 
processed  using  the  program  VuPoint  (S-Cubed  Inc.,  Division  of  Maxwell,  San  Diego,  CA).  Tlie  VuPoint 
program  produces  two  files — one  containing  calibration  infonnation  and  the  second  the  firing  data.  Under 
this  method  of  importing  data  to  BRLCB,  these  two  files,  which  arc  in  units  of  voltage  and  time,  are 
automatically  read  and  the  data  convened  to  pressure  data  utilizing  tlie  method  described  below. 

Conversion  to  Engineering  Units. 

(1)  To  convert  the  Nicolet  voltage  data  to  engineering  units,  the  first  .step  is  to  compute  the 
respoaso  factor  for  tlie  whole  clcctronic/digilal  data  acquisition  system  (i.c.,  tlie  amount  of  charge  [Q] 
represented  by  each  ADC  count),  'fliis  is  done  as  follows; 

(a)  Compute  charge  for  calibration  step 

Q  =  eVi  .  (3) 


where 


V-  =  calibration  voltage  input  to  the  charge  .'implifier* 

C  =  capacitance  of  cal  capacitor  in  charge  amplifier  {l.UUO  pF) 
Q  =  charge,  picocoulombs 

*Notc:  Tlie  V-  information  docs  NOT  appear  on  tlie  Nicolet. 
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(b)  Charge  per  ADC  count 


where 

F  =  response  factor  for  system  (picocoulombs/count) 

Q  -  charge  (picocoulombs) 

C'^  =  maximum  value  of  cal  step 
^mn  -  minimum  value  of  cal  step 

Note:  must  be  positive  and  negative. 

(2)  Converting  Data  to  Engineering  Units. 

(a)  Tlic  value  of  is  subtracted  from  each  pressure  count.  Tliis  amounts  to  a  baseline 
correction  to  get  a  corrected  ADC  count  for  the  pressure  at  each  point. 

(b)  Using  tlie  second-order  fit  equation  giving  pressure  in  tenns  of  charge  developed;  the 
pressure  is  computed  for  each  time  step,  as  indicated  below: 

P^A.B  (fC,,,)  *  C  (f  ,  (5) 


where 

P  -  pressure  (MPa) 

A,  B,  C  =  zero,  first-  and  second-degree  calibration  coefficients 

F  =  response  factor  (picocoulombs  per  ADC  count) 

=  counts  for  each  press  point 

2.5  Option  4.  Prepare  Finng  Information  File  (Program  MKINF.EXH).  In  Option  1 ,  Create  Master 
Infomiation  File,  inlbrmation  which  would  remain  constant  for  a  scries  of  experimental  firings  or 
simulations  was  entered  and  saved  in  a  file  termed  the  master  information  file.  In  this  option  (Option  4), 
that  information  specific  to  a  single  experimental  firing  or  simulation  is  apixmded  to  the  master 
infoimation  file  to  create  a  unique  file  known  as  the  .INF  file  or  information  file.  Since  the  infomiation 
file  is  created  from  a  master  infoimation  tile  created  in  Option  1,  die  first  query  posed  to  the  user,  after 
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Option  4  is  selected,  is  for  the  name  of  the  master  information  file,  as  shown  in  Tabic  16.  The  master 
information  file  name  should  be  given  with  both  the  drive  and  extension. 
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Table  16.  Request  for  Master  Infomiation  File 


Information  pertaining  to  an  individual  firing  is  entered  using  this  option. 
The  file  is  built  from  a  master  file  created  using  Option  1.  Enter  the  name 
of  the  master  file  which  will  be  used  in  creating  the  cunent  information 
file,  include  drive  and  extension. 


Next,  as  shown  in  Table  17,  the  type  of  computation  to  be  performed  is  determined.  Depending  on 
the  computation  mode  selected,  additional  infonnation  required  for  the  specific  computation  will  be 
requcittcd.  Specific  details  for  each  computational  mode  arc  provided  in  the  following  sections. 


Table  17.  Prompt  for  Selection  of  Computational  Mode 


Enter  the  type  of  computation  which  is  to  be  performed. 

1. 

Bum  rate  reduction 

2. 

Inverse  analysis  (generate  FTT) 

3. 

Surface  area  analysis 

4. 

Interrupted  burner 

5. 

ETC  reduction 

2.5.1  Bunt  Rate  Reduction.  To  perform  a  burning  rate  reduction  an  cxpcrimcnul  prcssurc-time  file 
is  required.  This  file  MUST  be  created  under  Option  3  beferc  Option  4  can  be  selected.  ITius,  when  the 
bum  rate  reduction  mode  is  selected,  tlic  first  request  is  for  the  name  of  the  pressure-time  file  created  in 
Option  3.  This  prompt  is  shown  in  Table  18.  After  the  pressure-time  file  is  examined,  tlic  user  is  asked 
to  verify  the  time  step. 


Table  18.  Prompt  for  Pressure-Time  File  Created  in  Option  3 


Enter  the  name  of  tlic  pressure-time  file  created  in  Option  3  which  is 
associated  with  tlic  information  file  being  created.  This  file  must  be 
ercated  before  the  information  file  is  created.  If  the  pressure-time  file  docs 
not  exist,  exit  the  program  and  complete  Option  3.  Remember  tlie  file 
name  is  given  without  extension  but  all  DOS  path  options  are  applicable. 


As  described  in  the  section  on  Option  3,  the  natne  of  the  pressure-time  file  (without  extension)  will 
determine  the  names  of  all  other  files  used  in  the  computation.  The  various  files  and  associated  names 
(tlie  name  consists  of  the  pressure-time  file  name  wiUi  a  specific  extension)  arc  shown  in  Table  19  as  they 
arc  p  esented  to  the  user  while  running  Option  4.  A  description  of  each  file  is  also  provided.  (Note:  In 
Table  19,  a.testl  is  the  example  test  name  used  in  Option  3,  this  portion  of  the  file  names  will  be  whatever 
name  was  used  for  the  pressure-time  file  created  in  Option  3.) 


Table  19.  Files  and  File  Names  Utilized  in  Bum  Rate  Reduction 


The  following  file  names 

will  be  used: 

Master  File 

a;testl.mas 

Created  in  Option  1 

Information  File 

a:testl.inf 

Created  now  in  Option  4 

Pressure-Time  File 

a;testl.pvi 

Created  in  Option  3 

Smoothed  Data  File 

a;tcstl.pdt 

Created  in  Option  5 

Output  File 

aitcstl.out 

Created  in  Option  6 

Graphics  File 

a;tcstt.dat 

Created  in  Option  6 

Bum  Rate  File 

a;tcstl.br 

Created  in  Option  6 

Pause  -  Please  enter  a  blank  line  (to  continue)  or  a  DOS  command. 

Master  File  (filcnamcl.mas):  File  created  in  Option  1  used  to  generate  the  information  file  currently 
being  created. 

Infonnation  File  (filename. inf);  File  currently  being  created. 

Pressure-Time  File  (filename .pvt);  Prcs.surc-timc  data  file  created  in  Option  3  associated  with  the 
information  file  being  created. 


Smoothed  Data  File  (filcnamc.pdi);  File  created  in  Option  5  which  will  contain  tlie  pressure-time  data, 
which  may  be  smoothed,  etc.,  directly  used  in  the  bum  rate  computation. 

Output  File  (filcname.out):  File  to  which  the  results  of  the  computation  will  be  printed. 

Graphics  File  (filcnamc.dat);  Special  file  created  to  use  with  graphic  packages  external  to  the  BRLCB 
program.  The  structure  of  this  file  for  all  the  computation  modes  is  shown  in  Table  20. 
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Table  20.  Smicturc  of  Graphics  File 


Computational  Option 

Column  1 

Column  2 

Column  3 

Bum  rate  reduction 
Inverse  analysis 

Surface  area 

Interrupted  burner 

ETC  reduction 

Pressure  (MPa) 

Time  (ms) 

Fraction  burned  (-) 
Pressure  (MPa) 
Pressure  (MPa) 

Bum  rate  (cm/s) 
Pressure  (MPa) 

Surface  area  ratio  (-) 
Bum  rate  (c/s) 

Bum  rate  (c/s) 

0.0 

0.0 

Depth  burned  (cm) 
0.0 

0.0 

Bum  Rate  File  (fileniune.br):  File  containing  results  of  the  computation  of  bum  rate  laws  in  the 
standard  r  =  bP"  form,  r  in  cm/s. 

After  the  file  names  arc  determined,  the  specific  information  concerning  the  hardware  and  charge  mass 
utilized  in  the  experiment  arc  requested.  The  necessary  information,  together  with  proper  units,  is 
presented  in  Table  21. 


Table  21.  Hardware  and  Charge  Information  Requested  in  Option  4 


Information 

Units 

Bomb  volume 

cm^ 

Initial  bomb  temperature 

K 

Propellant  mass 

g 

Igniter  mass 

g 

Propellant  temperature 

K 

Igniter  temperature 

K 

Finally,  information  concerning  the  pressure  range  for  the  bum  rate  laws  is  requested  (Table  22). 
Pressure  ranges  arc  based  upon  tlie  maximum  observed  pressure. 


'ITiis  completes  the  required  information  for  the  bum  rate  reduction  analysis,  the  infonnation  file  is 
then  written  and  the  program  returns  to  the  main  BRLCB  menu. 


Table  22.  Request  for  Pressure  Rjogcs  for  Bum  Rate  Laws 


One  ouqjut  of  the  bum  rate  analysis  is  a  scries  of  bP"  bum 
rate  laws.  The  program  has  the  option  to  determine  the 
bum  rate  laws  for  any  given  pres.sure  range  up  to  the 

Three  preset  ranges  can  be  used; 

5-10% 

10-25% 

25-75% 


1.  Use  predetermined  ranges. 

2.  Enter  new  ranges. _ 


2.5.2  Inverse  Analysis  (Generate  P/T).  Unlike  the  bum  rate  reduction,  the  inverse  analysis,  generation 
of  a  pressure  history,  docs  not  require  a  pressure-time  file;  thus.  Option  3  is  not  utilized  (Figure  2). 
Therefore,  upon  selecting  the  inverse  analysis  mode,  the  first  question  of  the  user  is  to  provide  the 
filename  to  be  used  for  the  infonnation  file  as  weU  as  all  the  other  files  created  during  the  analysis.  These 
files  and  their  names  arc  shown  in  Table  23,  Note  that  the  files  are  a  subset  of  those  used  in  the  burning 
rate  reduction  calculation  (Table  19).  (As  before,  a:testl  and  a:testpgen  arc  simply  example  file  names.) 


Table  23.  Files  and  File  Names  Utilized  in  the  Inverse  Analysis 


1  The  following  file  names  will  be  used: 

Master  File 

:  a;tcstl.mas 

Created  in  Option  1 

Information  File 

:  a;tcstpgen.inf 

Created  now  in  Option  4 

Output  File 

:  a;tcstpgcn.out 

Created  in  Option  6 

Graphics  File 

;  a:tcstpgen.dat 

Created  in  Option  6 

Pause  -  Please  enter  a  blank  line  (to  continue)  or  a  DOS  command. 

Next,  the  time  step  in  milliseconds  to  be  used  in  the  pressure  generation  is  requested.  As  in  the  burn 
rate  reduction,  the  user  is  then  prompted  for  specific  hardware  and  charge  information.  This  information 
is  identical  to  the  bum  rate  reduction,  sec  Table  21.  However,  unlike  the  bum  rate  reduction  calculation, 
to  generate  a  pressure  histoiy  the  propellant  linear  burning  rate  must  be  known.  BRLCB  allows  the 
burning  rate  to  be  entered  in  one  of  two  ways,  as  indicated  in  Table  24. 
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Tabic  24.  Methods  for  Entering  Propellant  Burning  Rate  Information 


Enter  the  method  by  which  the  bum  rate  information 
will  be  entered. 


1.  bP’'  bum  rate  lav;  one  law  for  each  layer 

2.  Table  of  pressure  vs.  rate 

Enter  your  choice  by  number  (1-2). 


In  Option  1  (ZjP"  bum  rate  law),  the  program  can  only  accommodate  one  law  per  grain  layer.  Under 
the  second  option,  tabular  pressure  vs.  rate,  no  such  restriction  is  imposed.  Bum  rates  at  intermediate 
pressures  are  obtained  through  linear  interpolation  of  the  table  of  log  (pressure)  and  log  (rate).  The 
information  file  is  then  written  and  the  program  returns  to  the  main  BRLCB  menu. 

2.5.3  Surface  Area  Analysis.  Required  input  for  the  surface  area  analysis  is  identical  to  the  bum  rate 
reduction  analysis,  except  that  information  concerning  the  propellant  burning  rate  is  required.  This 
information  is  obtained  as  in  the  inverse  analysis  (Table  24). 

2.5.4  Interrupted  But  ler.  Required  input  for  the  interrupted  burner  analysis  is  identical  to  the  bum 
rate  reduction  analysis. 

7.5.5  ETC  Reduction.  The  ETC  reduction  is  essentially  a  bum  rate  reduction  for  pressure-lime  data 
obtained  in  an  expenmental  firing  in  which  an  electrically  generated  plasma  has  been  injected  into  the 
closed  chamber  during  the  propellant  combusiion.  Thus,  in  addition  to  the  information  required  for  the 
bum  rate  reduction  analysis,  information  concerning  the  electrical  energy  injected  into  the  closed  chamber 
is  required.  BRLCB  expects  the  electrical  energy  information  to  be  contained  in  a  separate  file.  The 
prompt  for  this  file  name  (Table  25)  follows  after  the  pressure-time  data  file  is  examined  and  the  time  step 
verified. 


Table  25.  Prompt  for  Electrical  Energy  File 


For  the  ETC  reduction,  a  file  of  cumulative  electrical  energy  vs.  time  is 
required.  The  time  must  match  the  time  in  the  pressure-time  file,  the  units 
arc  seconds.  The  electrical  energy  is  in  units  of  megajoules  (MJ).  Enter 
the  file  name. 
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For  this  option  to  function  properly,  the  time  step  in  the  electrical  energy  file  must  be  the  same  as  the 
time  step  of  the  pressure-time  file  and  must  be  aligned  in  time  with  the  pressure-time  data.  Note  that  the 
file  is  to  contain  the  cumulative  electrical  energy  in  mcgajoules.  In  addition  to  the  electrical  energy  file 
name,  the  user  will  also  be  asked  for  the  total  electrical  energy  in  inegajoules  (MJ). 

2.6  Option  .5,  Smooth  Pressure-Time  Data  (Program  MKSMOOTH.EXE’).  In  Option  3  (Prepare 
Pressure -Time  Data),  pressure-time  data  arc  imported  into  the  BRLCB  program;  however,  no  attempt  is 
made  in  that  option  to  "massage"  or  prepare  the  data  for  the  desired  BRLCB  calculation.  Even  if  no 
modifications  to  the  pressure-time  data  imported  in  Option  3  is  desired  by  the  user,  Option  5  must  still 
be  executed  in  order  to  prepare  the  data  in  a  format  expected  by  the  data  analysis  program  (Option  6). 
Table  26  summarizes  the  data  massaging  available  and  indicates  whether  the  manipulation  will 
automatically  be  performed  or  is  optional. 


Table  26.  Data  Preparation  Options  Available  in  BRLCB 


Wildpoint  removal 

Optional 

Reduce  number  of  data  points 

Automatic 

Smoothing  data 

Optional 

Calculation  of  dP/dt 

Automatic 

Make  data  monotonic  increasing 

Optional 

The  option  begins  by  requesting  a  file  name.  As  discussed  in  Options  3  and  4.  this  file  is  the  file 
name  catered  in  Option  3  without  extension.  The  program  will  automatically  append  the  correct 
extensions  (.PVT  for  the  pressure-time  data  and  .INF  for  the  infonnation  file).  As  indicated  in  the 
previous  section,  this  option  will  create  a  file  with  the  same  file  name  but  with  an  extension  of  .PDl’. 
Once  the  file  name  is  obtained,  the  program  will  display  the  smoothing  and  differentiation  options,  as 
shown  in  Tabic  27. 
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Tabic  27.  BRLCB  Smoothing  and  Differentiation  Options 


Smoothing  and  Differentiation  Options 

1.  Use  a  fixed  bridge  length  (slivering  not  coasidered). 

2.  Use  a  "floating''  bridge  length  (slivering  not  considered). 

3.  Use  a  fixed  bridge  length  (slivering  considered). 

4.  Use  a  floating  bridge  length  (.slivering  considered). 

5.  Exit 

Please  enter  your  choice. 


Essentially,  smoothing  and  differentiating  is  performed  using  a  weighted  average  aty)ut  each  point. 
For  smoothing,  only  tlic  pressure  data  arc  utilized;  for  differentiation,  both  tlic  pressure  data  and  time  step 
arc  used.  The  number  of  points  used  in  the  averaging  is  called  the  "bridge  length."  In  BRLCB,  the  bridge 
length  can  vary  from  S  to  35  points,  but  mu.st  be  an  odd  number.  Coefficients  for  both  the  smootliing  and 
differentiation  are  listed  in  Appendix  M,  contained  in  volume  2  of  this  report.  To  illustrate,  consider  using 
a  5-point  bridge  Icngtli  for  smoothing  the  50th  prcs.surc  data  value,  P50.  Let  sl-s5  be  tlie  smoothing 
coefficients.  Ificn,  the  smoothed  pressure  value  which  will  rcpiac.,  the  50ili  pressua-;  point  is 


sl(P48)  +  s2(P49)  +  s3(P50)  +  s4(P51)  +  s5(P52)  . 


If  a  fixed  bridge  length  option  is  .selected  (Options  1  and  3),  then  Uiat  bridge  length  will  be  used  for 
all  the  data  except  f  or  data  points  at  the  beginning  and  end  of  the  data  record  where  the  bridge  length  is 
adjusted  so  that  indices  for  the  data  outside  the  data  bounds  arc  riot  requested.  For  example,  the  first  two 
and  la.st  two  data  points  are  never  smoothed;  data  point  3  is  smoothed  with  a  bridge  length  of  5;  data 
point  4  with  a  bridge  length  of  7;  etc.  Using  the  floating  bridge  length  options  (Options  2  and  4),  rcsulLs 
in  the  program  automatically  computing  the  bridge  length  to  u.sc  with  each  point.  This  bridge  Icngtli  will 
vary  with  each  point  and  is  selected  so  that  die  difference  between  the  highc.st  and  lowest  pressure  in  die 
smoothing  is  less  than  10%  of  the  maximum  ob.scrvcd  pressure  for  the  firing,  The  same  bridge  length 
computed  for  the  smoothing  is  uUlized  in  the  diffcrcritiauon.  Additional  details  on  the  .smoothing  and 
differentiation  methods  can  be  found  in  a  report  by  Doman  (1988). 

Slivering  refers  to  the  point  when  the  smallest  propellant  web  has  burned  through  and  is  often 
accompanied  by  a  distinct  change  in  the  slope  of  the  pressure-time  curve.  Smoothing,  using  data  points 
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on  both  sides  of  this  point  in  the  pressure-time  curve,  can  distort  the  data  which  will  be  reflected  in  the 
computed  burning  rates.  If  the  user  selects  either  Option  3  or  4.  .smoothing  and  differentiation  up  to  and 
after  tltc  slivering  pressure  are  performed  independently.  However,  BRLCB  docs  not  compute  the  pressure 
associated  witli  the  slivering  event.  Tltis  must  Le  supplied  by  the  user. 

(Note:  The  derivative  of  the  pressure,  dP/dt.  is  not  utilized  in  any  of  the  calculational  modes  of 
BRLCB.  However,  it  is  computed  since  this  information  is  utilized  by  some  researchers  and  propcUant 
formulators  to  determine  certain  propelhmt  characteristics,  namely,  "n'.lativc  quickness.") 

As  indicated  in  Table  26,  some  of  the  data  preparation  options  arc  automatic  while  others  are  optional. 
To  not  perform  cither  the  wildpoint  removal  or  smooiliitig,  enter  a  zero  when  the  program  requests  the 
number  of  wildpoint  or  smoothing  passes  to  perfonn.  (Note:  Multiple  wildpoint  removal  and  smoothing 
passes  arc  pemiitted  in  BRLCB.)  As  for  making  the  pressure  data  monotonic,  which  will  improve  run 
time,  make  tlie  desired  choice  when  prompted.  One  addition  option  is  available  for  the  wildpoint  removal. 
If  wildpoint  removal  is  selected,  the  user  is  able  to  indicate  a  tolerance  for  keeping  or  rejecting  a  point. 
Tlie  prompt  for  this  option  is  sluwn  in  Table  28. 

Table  28.  Prompt  for  Tolerance  on  Wildpr)lnt  Removal 

This  program  allows  the  user  to  enter  a  value  to  determine  what  tlic  cutoff 
will  be  for  discarding  or  keeping  wildpoints  or  outliers,  llie  value  entered 
can  be  any  number  greater  than  zero,  'flic  clo.scr  to  zero,  tlte  tighter  the 
tolerance  on  die  wildpoints. 

Lntcr  the  value  for  the  tolerance  on  Uie  wildpoints.  A  value  of  5  is 
generally  used. 


Essentially,  the  value  entcicd  can  be  thought  of  as  the  number  of  sttmdard  deviations  required  for  a 
point  to  be  considered  an  outlier  or  wildpoint.  Additional  details  can  be  found  in  Doman’s  report  (1988). 

nic  actual  flow  of  the  data  preparation  in  this  option  is: 

1.  Wildpoints  arc  removed  from  the  entire  data  set  if  tlie  option  is  selected.  The  number  of 
wildjx)int  passes  .selected  by  die  user  arc  perfonned  consecutively  at  this  point. 
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2.  The  number  of  data  points  are  reduced.  Points  are  deleted  from  the  beginning  of  the  file  so 
that  no  point  witii  a  pressure  below  80%  of  tiie  theoretical  igniter  pressure  is  retained.  All  points  after 
the  maximum  observed  pressure  arc  also  deleted.  A  summary  of  the  deleted  points  is  provided  to  the  user, 
see  Table  29. 


Table  29.  Summary  of  Deleted  Data  Points 


The  pressure  due  to  the  igniter  is: 

MPa 

The  starting  pressure  value  is: 

MPa 

N  points  have  been  deleted. 

This  corresponds  to  a  time  delay  of: 

ms 

The  indices  are: 

Pause  Please  enter  a  blank  line  (to  continue)  or  a  DOS  command 

Note:  If  the  ETC  bum  rate  computation  has  been  selected,  data  points  in  the  electrical 

energy  file  will  be  eliminated  in  an  identical  manner  to  maintain  tiic  time  alignment 
between  tlic  two  tiles. 


3.  The  reduced  data  file  is  smootli  with  the  number  of  passes  selected  by  the  user  according  to 
the  selected  smoothing  option. 

4.  The  pressure  derivative  is  computed  on  the  smoothed  data. 

5.  The  data  file  is  furtlicr  reduced  in  size  so  that  no  data  point  after  tlic  data  point  when  the 
theoretical  igniter  pres.surc  (with  heat  loss  considered)  is  first  obtained  has  a  pressure  below  tlic  igniter 
pressure.  (Note:  In  the  analysis,  the  igniter  is  considered  to  be  all  burnt  so  the  deleted  data  pxiints  would 
not  be  utilized  in  the  compulation.) 

6.  The  pre.ssurc  data  are  made  monoionically  increasing  if  that  opuon  ts  selected. 

7.  The  .PDT  created  by  this  option  is  written  with  the  format  given  in  Table  30.  However,  a 
maximum  of  999  points  will  be  saved  in  this  file.  If  points  are  deleted,  the  deletion  will  start  from  the 
beginning  of  the  data,  'fhe  reduction  in  points  is  by  UTincation,  not  decimation.  (Column  4  is  included 
for  the  ETC  reduction  only.) 
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Table  30.  Format  of  .PDT  File 


Column  1 

Column  2 

Column  3 

Column  4 

Time 

Pressure 

dP/dt 

Electrical  Energy 

(s) 

(MPa) 

(Million  MPa/s) 

(MJ) 

Note:  When  data  points  are  deleted,  die  time  is  adjusted  so  that  the  first  data  point  has  time  0.0  s. 


2.7  Option  6,  Peiform  Data  Analv.si.s  (Program  MK CAL. EXE).  As  the  name  indicates,  this  option 
performs  the  actual  data  analysis  indicated  in  Option  4,  Prepare  Firing  Information  File,  by  the  user.  The 
user  i.s  first  prompted  to  enter  the  file  name  (the  common  .INF,  .PVT,  and  .PDT  file  name)  associated  with 
the  analysis.  As  discussed  in  previous  sections,  the  file  extension  is  not  included.  (Note:  Entering  a  10 
at  the  file  name  prompt  will  allow  the  user  to  change  the  default  convergence  criteria  utilized  in  the 
program  for  computing  mass  burned.  The  current  default  is  lE-5.)  After  the  file  name  is  entered,  tlie 
program  will  access  aU  necessary  files  and  begin  the  analysis.  First,  the  program  will  determine  the  mass 
of  air  in  the  closed  chamber  based  upon  die  dau  contained  in  the  .INF  file.  The  user  may  elected  to 
ignore  the  air  mass  in  the  bomb  (c.g.,  the  closed  chamber  was  evacuated  before  the  firing)  and  is  provided 
Uiis  option,  as  shown  in  Table  31. 

Table  31.  Prompt  for  Air  Mass  Option 

The  computed  mass  of  air  in  the  bomb  is: _ g.  Do  you  wish  to  j 

change  this  value  to  0.0  g?  (yes  =  1,  no  =  2) 

Enter  your  choice. 


Next,  information  concerning  the  heat  loss  is  presented.  BRLCB  uses  the  term  heat  loss  fraction 
which  is  computed  as 

Heat  Loss  FracUon  .  1  -  ObscrvciJ  Maximum  Piessunt 

Theoretical  Maximum  Pressure 

If  the  analysis  mode  is  1  (Bum  Rate  Reduction),  3  (Suifacc  Area  Analysis),  or  5  (ETC  Reduction), 
the  heat  loss  information  will  be  presented  as  shown  in  Tabic  32.  In  general,  a  different  value  for  the  heat 
loss  fraction  should  not  be  entered.  The  information  is  provide  mostly  as  a  "sanity"  check.  Observed 
maximum  pressure  should  never  exceed  the  theoretical  maximum  pressure.  If  it  does,  the  input  values 
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to  the  program  should  be  carefully  examined.  On  average,  the  heat  loss  should  be  in  the  5%-10%  range 
(LOVA  propellants  up  to  20%),  heat  loss  fraction  of  .05-1.  When  heat  losses  greatly  exceed  these  values, 
the  input  data  to  the  program,  data  recording,  and/or  experimental  procedures  should  be  examined. 


Table  32.  Presentation  of  Heat  Loss  Infonnation 


♦■"Based  on  the  following  pressure  information; 
Observed  maximum  pressure  (MPa)  value  of  =  . 
Maximum  theoretical  pressure  (MPa)  value  of  = 
ITie  current  heat  loss  fraction  is:  _ 

Enter  a  different  value?  (yes  =  1,  no  =  2) 


However,  if  analysis  mode  2  (Pressure  Generation)  or  4  (Interrupted  Burner)  is  selected,  a  value  for 
the  heat  loss  fraction  must  be  entered.  The  prompt  to  enter  the  heat  loss  fraction  is  given  in  Table  33. 


Table  33,  Prompt  for  Heat  Loss  Fraction,  Analysis  Modes  2  and  4 


For  the  pressure  generation  option  (Option  2)  or  the  interrupted  burner 
option  (Option  4),  a  heat  loss  fraction  must  be  entered.  The  heat  loss 
fraction  is  in  the  range  of  0.0-1. 0.  For  example,  a  value  of  0.1  means 
10%  of  the  total  energy  will  be  considered  lost  as  heat  to  Ihc  chamber  wall 
during  the  calculation. 

Enter  the  value  for  the  heat  loss  factor.  A  decimal  between  0.0  and  1.0. 


After  air  mass  and  heat  loss  arc  handled,  tlic  analysis  is  started.  As  mentioned  earlier,  the  model 
assumes  that  the  igniter  is  all  burnt  at  tlic  start  of  the  analysis.  Thus,  points  at  tlie  beginning  of  the  data 
flic  may  be  ignored.  Information  concerning  the  ignored  data  points  is  provided  to  the  u.scr,  as  shown  in 
Table  34. 
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Table  34.  Information  Concerning  Data  Points  Ignored  in  tlie  Calculation 


For  the  calculation  to  be  performed  correctly,  all  pressures  must  be  above 
the  igniter  pressure.  For  the  current  data  set  points  had  to  be  deleted  to 
obtain  all  pressures  above  the  igniter  pressure.  Information  relative  to  the 
deleted  data  points  is: 

Number  of  points  deleted: 

Time  interval  of  deleted  points  (ms): 

Igniter  pressure  (MPa): 


2.8  Option  7.  Prepare  Output  (I^rogram  MKOtJT.EXE).  This  option,  as  tltc  name  indicates,  is  used 
to  prepare  both  hard  copy  and/or  magnetic  output  for  the  analysis.  The  option  provides  for  computation 
of  bum  rate  laws  and  graphical  output  as  well  as  a  one  page  summary  sheet  for  bum  rate  reductions 
(analysis  modes  1 , 4,  and  5). 

Appendix  C  contains  a  sample  output  Irom  BRLClt.  First,  the  input  infomiaiiori  provided  by  the  user 
is  echoed.  Next  is  a  complete  listing  of  all  data  u.sed  in  the  computation  and  information  concerning 
maximum  chamber  properties,  together  with  heat  loss  information.  Tabular  data  from  die  computation 
follow.  Table  35  summarizes  the  infonnation  for  the  nine  columns  of  the  tabular  data. 


Table  35.  Structure  of  Tabular  Data  Output  for  BRLCB 


Column 

Infonnation 

Units 

1 

Computational  step 

_ 

2 

Propellant  layer 

- 

3 

Time 

ms 

4 

Chamtx'-r  pressure 

MPa 

5 

Unbumed  propellant  mass 

g 

6 

Bum  rate 

enVs 

7 

Propellant  surface  area 

cm^ 

8 

Average  chamber  temperature 

K 

9 

Depth  burned 

cm 

This  option  also  allow.s  the  user  to  perform  a  Fast  Fourier  Transform  (FFI’)  analysis  of  computed 
burning  rates  if  the  rates  display  a  large  amount  of  oscillations.  Trial  and  error  will  be  required  to 
determine  tlic  best  cut-off  value  for  the  m’  frequency.  Generally,  selecting  a  low  frequency  will  provide 


the  best  results.  Figure  19  shows  burning  rates  computed  for  an  M5  propellant.  As  can  been  seen  in  the 
figure,  some  oscillation  of  the  burning  rate  is  present.  After  performing  a  low  pa.ss  filter  on  the  data,  the 
burning  rate  shown  in  Figure  20  is  obtained  which  displays  much  less  oscillation.  A  comparison  of  the 
two  burning  rates  is  provided  in  Figure  21.  In  the  program,  the  user  is  able  to  determine  the  effect  of  the 
ITT  and  decide  whether  to  use  the  original  or  filter  burning  rate  data. 


2.9  Sample  Problem.  On  the  diskette  containing  the  source  files  for  BRLCB  is  a  subdirectory  labeled 
SAMPLE.  This  directory  contains  all  the  files  from  a  sample  BRLCB  burning  rate  reduction  (see 
Table  36).  Tfie  information  can  be  used  either  to  verify  the  computational  portion  of  the  BRLCB  program 
or  provide  the  u.scr  with  a  ''practice"  set  of  data  to  utilize  in  learning  how  to  run  the  entire  BRLCB 
program  for  burning  rate  reduction. 


Table  36.  Files  in  Subdirectory  SAMPLE 


testi.pt 

TESTl.MAS 

TESTl.PVT 

TESTl.INF 

testi.br 

TEST1.PDT 
TEST  LOUT 
TEST1.DAT 
TESTl.PRT 


To  verify  the  computational  portion  ol'  tlie  program,  follow  die  following  procedure. 

1.  Copy  die  files  TESTl.INF,  TESTl.PDT,  and  TESTl.BR  to  a  di.sk  in  drive  A  (diesc  files  mu.st 
be  on  drive  A). 

2.  Change  directory  to  the  directory  containing  die  BRLCB  program. 

3.  Start  die  program  by  typing  BRL  and  select  Option  6. 

4.  At  die  file  name  prompt,  enter  AnTiSTl,  accept  die  air  mass  as  computed,  and  set  the  heat 
loss  to  0.0. 
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Figure  19.  Original  unfiltered  M5  burning  rate  data. 
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Figure  2 


5.  When  the  program  returns  to  the  main  menu,  select  Option  7.  Respond  to  Uie  file  name 
prompt  again  with  A:TEST1  and  do  not  perform  any  filtering  on  the  bum  rate  data. 

The  printed  output  can  be  compared  with  the  information  in  the  file  TESTl.PRT.  If  these  files  are 
not  identical,  repeat  the  process. 

Instead  of  simply  verifying  the  computational  portion  of  the  program,  it  is  recommended  that  the  user 
exercise  all  the  options  of  the  program  by  using  the  sample  data  to  perform  the  burning  rate  reduction 
from  the  beginning.  In  this  case,  only  one  file  (TESTl.PT)  needs  to  be  copied  to  drive  A  (or  to  the 
directory  where  the  BRLCB  executable  programs  arc  located).  Then,  the  user  should  mn  all  options  (1-7) 
except  Option  2,  in  order.  Necessary  information  for  the  reduction  is  listed  in  Table  37.  Results  can  be 
compared  with  the  file  TESTl.PRT. 


Table  37.  Information  Required  for  Sample  Burning  Rate  Reduction 


Propellaiit: 

Grain  geometry 

Sphere 

Diameter 

0.072  cm 

Impetus 

1,165.4  J/g 

Flame  temperature 

3,680  K 

Molecular  weight 

26.25323 

Gamma 

1.2155 

Covolume 

0.966  cm^/9 

Density 

1.6  g/cm^ 

Mass 

13.1  g 

Initial  temperature 

294  K 

Igniter: 

Same  Uicnnochemistry  as  propellant 
Mass 

0.5  g 

Initial  temperature 

294  K 

Hardware: 

Bomb  volume: 

65.5  cm^ 

Time  Step: 

0.025  ms 
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3.  THEORETICAL  ANALYSIS  OF  CLOSED-CHAMBER  COMBUSTION 


A  closed  chamber  is  a  rather  simple  device  which  monitors  transient  combustion  with  a  single 
diagnostic — a  wall-mounted  pressure  transducer.  Except  for  the  initial  conditions,  the  pressuie-timc  record 
is  the  only  infonnation  available  to  deduce  an  effective  linear  regression  rate  of  the  material  confined  in 
the  chamber.  Of  course,  predicting  this  regression  (or  burning)  rate  is  not  a  new  problem.  A  number  of 
methods  of  varying  levels  of  sophistication  have  been  advanced  over  the  years,  as  illustrated  in  Robbins 
and  Horst  (1976);  Price  and  Juhasz  (1977);  Oberle,  Juhasz,  and  Griffic  (1987);  Robbins  and  Lynn  (1988); 
and  references  therein.  Although  a  comprenehsivc  review  is  beyond  the  scope  of  the  present  report,  it 
appears  that  two  basic  approaches  have  been  employed.  The  first  method  derives  a  closed-form  solution 
for  mass  remaining  (or  consumed)  in  the  chamber  as  a  function  of  the  experimental  value  of  pressure, 
The  second  method  derives  a  differential  equation  for  tlic  rate  of  change  of  propellant  mass  which  must 
be  numerically  integrated  in  time.  This  differential  equation,  of  course,  must  be  accompanied  by  some 
estimate  of  the  time  derivative  of  chamber  pressure  from  the  experimental  data. 

In  the  authors’  opinion,  tlie  first  mctliod  (closed-form  solution)  is  superior.  The  second  method  can 
be  encumbered  by  stability  and  accuracy  problems  associated  with  forward  time  integration  and  also 
requires  that  the  experimental  pressure-time  data  be  differentiated.  With  an  objective  to  eliminate  as  many 
sources  of  error  as  possible,  the  present  analysis  will  be  based  on  a  closed-form  solution  (first  method). 
Although  derived  independently,  the  analysis  shares  the  general  approach  first  outlined  by  Robbins  and 
Horst  (1976). 

The  closed  chamber  analysis  developed  here  invokes  several  assumptions  common  to  a  "well-stirrcd 
reactor."  Velocities  within  tlie  chamber  arc  assumed  vanishingly  small,  and,  hence,  balance  of  momentum 
implies  spatially  uniform  pressure.  For  consistency,  kinetic  energy  is  assumed  negligible  compared  to 
stored  thcima!  and  chemical  energy.  All  otlicr  properties  of  the  gas  phase  within  the  chamber  are  spatially 
invariant  as  the  result  of  tlic  "well-mixed"  assumption.  However,  solid-phase  properties,  as  discussed 
below,  can  vary  as  a  function  of  1)  the  particular  segment  of  propellant  which  is  burning  and  2)  how  much 
mass  of  that  segment  has  been  consumed.  Thus,  properties  of  the  combustion  gases  entering  the  chamber 
at  any  given  time  will  depend  upon  the  amount  of  solid  material  which  has  been  consumed  by  tlic 
combustion  process,  llic  chamber  initial  condition  assumes  that  the  igniter  material  has  been  consumed, 
although  a  description  of  igniter  combu.stion  could  easily  be  added  in  the  future. 
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3.1  Variable  Property  Description.  The  present  closed-chamber  combustion  analysis  is  not  restricted 
to  homogeneous  solid  propellant  material  with  spatially  uniform  properties.  Provisions  are  made  for 
propellant  properties  such  as  density,  stored  thermal  energy,  etc.  to  have  a  specified  variation  as  a  function 
of  the  depth  beneath  the  initial  surface.,  which  is  related  to  the  amount  of  material  which  has  burned  away. 
Although  the  propellant  properties  must  be  time  invariant,  the  description  of  spatial  variability  has 
substantial  flexibility  while  at  the  same  time  incorporating  the  previous  description  of  layered  propellant 
(Oberlc  and  Kooker  1989).  Consider  the  schematics  shown  in  Figures  22-24,  where  v  represents  any 
propellant  property  and  x  represents  a  coordinate  which  is  zero  at  the  initial  surface  and  taken  as  positive 
in  the  direction  of  the  interior  of  the  propellant  grain.  The  propellant  material  is  then  subdivided  into  N 
layers  or  segments  in  the  x  coordinate,  where  the  segments  do  not  have  to  have  equal  thickness.  Let  a 
subscript  1  denote  the  value  associated  with  the  beginning  of  the  segment  and  a  subscript  2  denote  the 
value  at  the  end  of  the  segment.  A  basic  assumption  of  the  analysis  is  that  property  variations  wiu  hi  a 
given  segment  arc  linear  in  the  x  coordinate.  Between  segments,  however,  no  continuity  is  assumed  or 
required;  both  the  function  and  its  spatial  derivatives  may  be  discontinuous.  The  flexibility  provided  by 
this  description  is  illustrated  in  Figures  22-24.  Tlie  continuous  property  variation  shown  in  Figure  22  is 
enforced  by  setting 


¥i  =  ¥i  for  I  <  i  <  N  -  I  , 

(8) 

where  the  accuracy  of  the  representation  by  linear  segments  wiU  be  governed  by  the  chosen  thickness  of 
the  .segments.  A  material  with  distinct  layers  (see  Figure  23),  each  of  whicli  has  uniform  properties,  is 
described  with 


¥i  =  ¥9  for  1  ^  i  N  . 

And,  as  sketched  in  Figure  24,  the  analysis  will  also  describe  a  material  composed  of  a  combination  of 
distinct  layers,  each  of  which  has  variable  properties. 

An  analysis  of  closed-chamber  combustion  must  relate  chamber  pressure,  at  any  instant  of  time,  to 
the  amount  of  material  consumed  by  the  combustion  process.  This  requires  some  type  of  "progress 
variable"  to  monitor  the  extent  of  completion  of  the  combustion  process.  An  additional  role  for  this 
variable  is  to  link  the  in.stantancous  value  of  burning  surface  area  to  the  amount  of  solid  propellant 
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Figure  22.  Continuous  prope.rty  variation. 
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mass  which  has  been  consumed.  The  present  analysis  will  choose  the  progress  variable  to  be  the  fraction 
of  propellant  mass  in  segment  "i"  remaining  in  the  chamber,  denoted  as  4,-  and  defined 


4.- 


(10) 


where 

=  solid  propella»^t  mass  initially  contained  in  segment  "i" 

=  solid  propellant  mass  of  segment remaining  in  chamber. 


Note  that  when 

the  segment  is  unbumed  ==>  =  1 

the  segment  is  consumed  ==>  =  0. 

Clearly  then,  the  fraction  of  propellant  mass  of  segment  "i"  which  has  been  consumed  is 

1  - 

'I'his  (juantity  is  directly  related  to  the  depth  of  material  which  has  been  consumed  by  the  combustion 
process;  tlic  relationship  is  given  by  an  integral.  Assume  that  a  differential  volume  element  can  be 
represented  as 


dV  =  A(x)  dx 

where  "a:"  is  the  depth  coordinate  discussed  above  and  A(x)  is  tlie  grain  surface  area  at  depth  ''x".  Note 

that  the  direction  of  coordinate  x  is  always  normal  to  the  surface  A(x).  Now  let  x,  denote  the  depth  at 

‘i 

which  segment  "i"  begins,  and  X2.  denote  the  depth  at  which  it  ends.  When  combustion  of  segment 
is  underway  (i.c.,  Xi,<  \  <  X2),  then  the  mass  of  segment  "i"  which  has  been  consumed  is  simply 

X 

j  p/z)A(z)dz  . 

-*1. 

I 


■ 
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Then  the  relationship  tjctween  1  -  (the  fraction  of  propellant  mass  in  segment  "j"  which  has  been 
consumed)  and  x  -  x,  (the  depth  burned  of  segment  "i")  is  defined  by  the  mass  balance, 

JC 

(1  -  -  J ()^(z)A(z)d2  .  (11) 

Xx 


This  represents  the  crucial  link  between  grain  geometry  (depth  burned  and  surface  area),  variable 
propellant  density,  and  mass  consumed  by  the  combustion  process.  In  practice,  given  a  value  of  an 
iteration  procedure  is  used  to  determine  the  corresponding  value  of  "j:"  which  balances  Equation  (11)  to 
within  some  tolerance.  Also  note  the  obvious  equality  when  x  =  ,  i.c.. 


m 


"  -  J  p^(z)A{z)dz  . 


Now  the  assumed  linear  dependence  of  any  property  'E  can  be  expressed  as 


Ti.  ^  (4^2.  -  Ti.).iil/l>\. 
‘  ‘  ‘  (^2.  ■  ^1.) 


(12) 


(13) 


Because  the  present  analysis  must  allow  for  solid  propellant  with  variable  properties,  the  balance  equations 
for  combustion  will  be  a  function  of  certain  average  values  which  account  for  the  "history"  of  the  property 
variation  to  tltat  lime,  The  average  value  of  a  quantity  (in  the  i'*’  segment)  which  ha.s  been  consumed  i.s 
simply 

X 

j  H',  (7)^(z)dz 

^■(x)  =  :!ii_ -  (14) 

j*  A{z)dz 
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The  average  value  of  a  quantity  (in  the  i*  segment)  remaining  in  the  chamber  (not  vet  consumed)  is  then 


^2. 

J  4',  (2)A(z)«iz 

-  (15) 

•*2|- 

J*  A(z)  dz 

X 


Note  that  if  the  quantity  'V{x)  is  uniform  (constant  value,  e.g.,  '?•)  in  the  "ith"  segment,  as  shown  in 
Figure  23,  then  Equations  (14)  and  (15)  reduce  to 


(16) 

Thus,  in  the  special  case  of  constant  property  layers.  Equation  (16)  defines  both  the  average  value  of  the 
quantity  which  has  been  consumed  and  the  average  value  of  the  quantity  remaining  in  the  chamber. 

3.2  Governing  Equations.  The  equations  governing  combustion  of  solid  propellant  in  a  closed 
chamber  arc  straightforward.  Since  no  mass  enters  or  leaves  the  chamber,  conservation  requires  that  any 
solid  mass  consumed  in  the  combustion  process  must  reappear  in  the  gas  phase.  Since  the  velocities  of 
botli  phases  aie  assumed  negligible,  a  momentum  balance  in  the  chamber  reduces  to  a  statement  of 
spatially  uniform  pressure.  Tlic  balance  of  energy  monitors  the  transfer  of  stored  chemical  energy  in  the 
solid  to  the  gas  during  combustion,  and  accounts  for  any  heat  loss  through  the  chamber  walls.  Finally, 
a  gas-phase  equation  of  state  must  relate  pressure,  temperature,  and  gas  volume. 

Some  notation  is  required  to  express  these  balance  laws.  For  a  chamber  of  fixed  volume,  let 

P^f^  s  chamber  pressure, 
s  chamber  temperature, 

Kh  s  closed  chamber  volume 

^fres  =  chamber  volume  not  occupied  by  solid  material. 
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Then  define 


b  m  covolumc, 

C^,  m  specific  heat  at  constant  volume,  [here,  S?  /  (7  -  1)] 
e  3  specific  internal  energy, 
m  »  solid  propellant  mass  remaining  in  chamber, 

M  s  molecular  weight, 

(2^,  3  cumulative  heat  loss  to  chamber  walls 

E  universal  gas  consuint/molecular  weight  =  R^j  M, 

y 

p  E  density. 

Define  the  subscripts 

(  =>  pertains  to  solid  phase  (solid  propellant), 

(  =>  pertains  to  gas  phase  (combustion  products), 

(  ),^=>  pertaias  to  gas-phase  igniter  products, 

(  )„  =>  {xtrtains  to  air  initially  in  chamber, 

(  =>  pertains  to  chamber  conditions, 

(  )j  =>  evaluated  at  the  beginning  of  the  segment, 

(  =>  evaluated  at  the  end  of  the  segment, 

(  ),  =>  pertains  to  the  itli  segment  or  layer  of  the  solid  propellant. 

A  superscript  "o"  means 

(  )°  =>  pertains  to  the  initial  condition  (at  time  zero). 


For  each  segment  i,  the  initial  conditions  provide  rn^ ,  the  initial  mass  of  propellant  in  segment  i.  At 

i 

any  otlier  time,  is  the  mass  of  propellant  in  segment  i  remaining  in  the  chamber  as  solid.  Tlic  ratio 
of  these  two  quantities  is  chosen  as  the  progress  variable  Thus  by  definition, 
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'"iXO  »  •  ^,(0 


(17) 


and  mass  conservation  for  each  propellant  segment  requires  that 

m^XO  -  m  "  -  m^XO  " 'm/  [l  “  ^,(0]  ■  (18) 

Subscript  i  denotes  properties  associated  with  one  of  tlic  N  segments  of  the  solid  propellant  material.  It 
is  assumed  that  the  segments  are  numbered  sequentially,  with  the  segment  labelled  i  =  1  at  the  outer 
surface  of  the  unbumed  grain.  Note  that  at  any  time  before  segment  j  begins  to  bum,  is  unity:  at  any 
time  after  segment  i  has  been  consumed,  is  zero. 

Expres.sions  for  the  gas-phase  equation  of  state  and  the  energy  balance  f  or  the  chamber  will  dejxmd 
upon  certain  time-varying  average  quantities;  i.c.,  thermochemical  and  physical  properties  of  tlic 
proijcllant.  nicse  average  values  at  time  t  are  functions  of  tlie  solution  at  time  f.  Clearly,  the  desired 
solution  will  require  an  iteration  procedure.  Denne  tlic  following  average  quantities  which  are  functions 
of  the  (inal  solution  at  time  /: 


Pi'.  ”  "‘j.  I  (Average  density  of  material  not  yet  consumed.)  , 


which  employ  Equation  (15),  and 


(19) 


=  m°  <■  (Average  specific  internal  energy  of  material  not  yet  consumed)  ,  (20) 

i  i 


b  =  *  (Average  covolumc  of  material  consumed.)  , 


(21) 


C^,  =  *  (Average  specific  heat  at  constant  volume  of  material  consumed.)  ,  (22) 


(Average  value  of  universal  gas  constant/molccular 
|_  weight  of  material  consumed.) 
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(23) 


■ratWiCTBtttTWWlliS^^ 


which  employ  Equation  (14).  It  is  important  to  note  that  Equations  (19)  and  (20)  define  tlic  average  of 
a  quantity  not  yet  consumed,  whereas  Equations  (21)-(23)  define  the  average  of  a  quantity  which  has  been 
consumed.  The  reader  should  note  that,  as  defined  above,  the  "over  bar"  quantities  arc  all  functioas  of 


A  covolume  equation  of  state  for  the  gas  mixture  can  be  written, 


^ch  \yfree  ~  voluHic  occupicd  by  gas  molcculcsj 
N 

-  T-h  ^  nmoles  of  gas)  . 

“  i  -  1 


(21) 


The  volume  of  the  chamber  not  occupied  by  the  solid  propellant  (Vy^ee)  simply 


N 

eh-  ^  P^.  ■  4/  • 
1-1  ‘ 


(22) 


Then  the  equation  of  .state  [Equation  (21)]  becomes 


eh 


N 


^ch-  ^  P».  ~  2:  bg.(l-^,.) 

i  -  1  ‘  i  =  1  ' 


=  T,. 


eh 


N 


9?;^  (1  -  ^,) 

i  -  1  ' 


(23) 


For  convenience,  define  tlie  tcims 


N 


=  '^Ch  -  Z  b  , 


i  =  1 


^2  =  ^ig  ^ 

J  =  I  ' 
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neither  of  which  arc  constant;  both  must  be  updated  in  the  iteration  procedure.  Equation  (26)  can  then 
be  written, 


N 

M  1 

p 

^ch 

/Cl  -  L 

J  =  1 

(p^.  - 

-Teh' 

(27) 


To  write  the  energy  balance  in  the  chamber,  let  E  represent  the  total  energy  contained  in  Then 
the  first  law  of  thermodynamics  states  that 


£  ^  -  £(r)  +  Q^U)  (28) 

where  QJt)  is  the  cumulative  heat  loss  from  the  chamber.  For  tlic  ETC  reduction,  Qjt)  =  cumulative 
heat  low-cumulative  electrical  energy  input.  The  initial  energy  is  simply. 


£" 


u 


N 


(29) 


which  neglects  stored  thermal  energy  in  the  solid  material  at  tlic  initial  temperature  compared  to  stored 
chemical  energy.  For  any  time  greater  than  zero,  the  igniter  material  has  been  conveited  to  gaseous 
products,  and  the  total  energy  can  be  represented  as 


N 

E(t)  =  E 
i 


—  0 
=  1  ‘ 


eh 


a 


N 

^  C,  (1  -  4,) 
i  =  1 


(30) 


Hence,  the  first  law  in  Equation  (28)  becomes 


E°  -  Q^O) 


N 
S 
=  1 


—  0 
I 


'  eh 


^  .(1-4,) 

'«  “  i  -  \ 


(31) 


Again  for  convenience,  define  the  tenn 

ATj  £  m- 

'■«  /  =  1  Xi 
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which  also  must  be  updated  in  the  iteration  procedure.  Solving  Equation  (31)  for 


T.. 


eh 


-QJt)  - 


N 

K^-  T. 
i  «  1 


(32) 


Tlic  analysis  which  follows  is  based  on  tlic  convention  that  propellant  segments  arc  labeled  i  =  7,...,N 
in  the  order  in  wliich  combustion  will  reach  them.  Thus,  when  segment  k  is  burning  U<k<!^< 

a)  segments  7 (k-1)  have  been  consumed  «  0],  and  b)  segments  (k+1) N  are  unbumed 

[^'*5  »  1 1 .  Witli  this  convention  in  mind,  define  the  following  terms: 

N 

i  *  k 

N 

V'  . 

i  ^  k 

N 

5i  . 


C^^E° 


2.(0 


N 
=  1 


Note  that  C;,  C2,  Cj,  and  Q  must  also  be  updated  in  the  iteration  procedure.  The  gas  mixture  equation 
of  state,  Equation  (27),  can  then  be  written  as 


(c,  -  (p.^  -  ^,,1  ^4  -  T.,  {c,  -  5,^  54  , 


(33) 
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and  the  chamber  energy  equation,  Equation  (32),  becomes 


T 


eh 


c,  -  4, 

Sk 


(34) 


Substituting  Equation  (34)  into  Equation  (33),  assuming  that  is  given  by  experimental  data  from 

tlic  chamber,  the  result  is  a  quadratic  equation  for  the  fraction  of  propellant  mass  in  segment  k  which 
has  yet  to  be  consumed. 

Qa^I  +  QB^k  +  Gc  =  0 


where 


—  O 


Gfl  “  Pch 


Cl 


Sk 


Qc  ®  C2C4  -  ^7,  • 


The  solution  of  Equation  (35)  is  simply 

-  -  {c!o  +  (qb  -  ‘I  e^Cc }  /2  s.,  . 

This  value  of  is  used  to  update  all  the  time-varying  average  quantities  as  well  as  the  terms  (e.g.,  Kj, 
K2,  Kj)  which  are  fimctions  of  these  quantities.  Then  a  new  solution  for  is  obtained  from 
Equation  (36),  and  tlie  process  is  repeated  until  successive  iterates  are  equal  to  within  some  tolerance. 
Experience  with  the  test  eases  has  shown  convergence  to  be  quite  rapid. 


Once  the  solution  for  the  has  been  determined,  chamber  tempxirature  follows  directly  from 
Equation  (34).  Note  that  the  value  of  chamber  pressure  at  which  the  1^^  segment  "bums  out"  is  given 
explicitly  hy  Qq  =  Q  or 
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(37) 


p  I  , 

I  burnout  ^ 

Once  the  segment  is  consumed  (fj^  =  0  and  hence  =  0).  the  analysis  automatically  assumes  that 
the  <:  +  7  segment  begins  burning  with  the  initial  condition  ^  ^  or  +  i  ®  1 .  When 

the  last  or  AT*  segment  has  been  coasumed,  the  function  m/t),  the  total  propellant  mass  remaining  in  the 
chamber, 

N 

m,{t)  -  E  m/  ^,(0  .  (38) 

i  =  1  ‘ 


has  been  determined.  The  time  derivative  dm^ldt  can  easily  be  detennined  by  finite  difference  with 
second-order  accuracy.  In  several  checkout  cases,  this  finite  difference  derivative  was  compared  to  the 
differential  of  a  cubic  spline  and  found  to  have  the  same  accuracy.  Given  the  instantaneous  surface  area, 
A(t),  as  a  function  of  depth  burned,  the  effective  linear  regression  rate  is  simply 


dm^ldt 

"mTo 


(39) 


It  is  well  known  that  the  values  of  linear  burning  rate  deduced  from  a  closed  chamber  pressure  record 
will  be  influenced  by  the  amount  of  heat  (or  energy)  lost  to  the  chamber  boundaries  or  walls.  In  the 
present  analysis,  this  is  represented  by  the  term  Q^,  the  cumulative  heat  loss  to  the  chamber  walls. 
Possibly  because  Uic  computation  of  is  not  trivial,  there  seems  to  be  no  general  agreement  on  the 
calculation  procedure.  In  many  cases,  however,  the  calculation  depends  upon  Q,„^  the  amount  of  heat 
lost  to  the  chamber  walls  up  to  the  time  of  maximum  pressure  [currently  in  BRLCB,  is  computed  as 
/  ^cbnuix)  *  Qmaxi-  Qmax  con^puted  in  the  following  way;  Equation  (26)  gives  a  value  of  7^.^  based 
on  the  measured  value  of  maximum  pressure;  Equation  (31)  is  then  solved  for  the  coniesponding  value 
of  heat  loss, 
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4.  VALIDATION 


To  validate  thie  computer  implementation  of  the  mathematical  model  and  the  various  grain  form 
functions,  a  series  of  test  cases  were  performed.  The  test  cases  can  be  divided  into  five  scries  with  the 
following  objectives. 

Scries  I.  Validate  the  computer  implementation  of  the  mathematical  model  'Jirough  comparison  with 
analytically  generated  data. 

Scries  II.  Quantify  the  affect  of  reduced  accuracy  in  the  input  pressure  data  on  the  accuracy  of  the 
deduced  burning  rate. 

Scries  III.  Validate  the  grain  form  functions. 

Scries  IV.  Validate  the  ETC  Burning  Rate  Reduction  option. 

Scries  V,  Validate  the  computer  implementation  for  multilayered  propellants  witli  comslant  properties 
and  varying  properties. 

4. 1  Scries  1  Test  Cases.  To  validate  the  compuUT  implementation  of  the  mathematical  model,  results 
from  BRLCB  arc  compared  to  an  analytic  (no  numerical  methods  utilized)  solution  of  the  closcd-chimbcr 
problem  (valid  only  for  single  perforated  grain  with  a  bum  rate  exponent  of  i)  using  a  procedure  provided 
by  Robbins  and  Lynn  (1988).  Table  38  summarizes  the  input  parameters  used  to  generate  the  analytic 
data  as  well  as  the  inputs  to  the  various  options  of  BRLCB.  Sufficient  information  is  available  from  the 
analytically  generated  data  to  validate  tlic  burning  rale  reduction  option  (1)  [and,  hence,  the  iniermptcd 
burner  option  (4)],  the  pressure  generation  option  (2),  and  the  surface  area  analysis  option  (3).  The  ETC 
bum  rate  reduction  is  validated  in  test  scries  IV. 

Tc.st  Case  1.  In  this  tc.st  case,  the  pressure-lime  data  from  the  analytic  solution  of  Robbins  and  Lynn 
arc  used  to  validate  the  burning  rate  reduction  (and  thus  the  inicmjptcd  burner  option)  option  of  BRLCB. 
RcsulLs  from  the  BRLCB  bum  rate  analysis  are  given  in  Figures  25  and  26.  The  deduced  burning  rate 


Table  38.  Inpul  Parameters  for  Robbins  and  Lynn  Analytic  Solution 


Qosed-chamber  volume: 

205.002  cm^ 

Grain  geometry: 

Single  Perforation 
length:  2.54  cm 
diameter:  .635  cm 
perf  diam:  .127  cm 

Propellant: 

Impetus: 

1,009.5959  J/g 

Flame  temperature: 

3,400  K 

Ratio  of  specific  heats: 

1.2299 

Molecular  weight: 

28 

Covolume: 

.90318  cm^/g 

Density: 

1.6608  g/cnr 

Total  mass: 

56.69905  g  +  .09072  g  igniter 

Bum  rate  law: 

r  =  .73679  P^  °  cm/s 

is  shown  in  Figure  25  with  the  corresponding  percent  error  with  the  assumed  burning  rate  given  in 
Figure  26.  As  can  be  seen  from  the  figures,  the  deduced  burning  rate  is  linear  (on  log-log  scale)  with  a 
maximum  percent  error  of  approximately  0.002%.  Thus,  it  appears  that  the  computer  implementation  is 
accurately  predicting  propellant  bum  rates  from  closed-chamber  pressure-time  data. 

Test  Case  2.  Using  the  same  inputs  (Table  38)  as  in  the  Robbins  and  Lynn  model,  BRIXB  was  used 
to  generate  the  closed-chamber  pressure-time  trace  shown  in  Figure  27.  A  comparison  with  the  Robbins 
and  Lynn  generated  pressure  history  (in  terms  of  percent  error)  is  provided  in  Figure  28.  Although  the 
maximum  percent  error  has  increased  by  an  order  of  magnitude  (0.02%  vs.  0.002%)  compared  to  test 
case  1,  the  authors  believe  that  this  is  still  within  the  round-offAruncation  error  associated  with  the 
numeric  methods  employed  in  BRLCB. 


Test  Case  3.  As  in  test  ease  1,  the  analytic  pressure  history  is  used  as  input  to  the  BRLCB  code. 
However,  in  this  test  ease,  the  assumed  bum  rate  (Table  38)  is  used  to  determine  the  grain  surface  area 
profile.  Thus,  the  surface  area  analysis  option  of  BRLCB  is  being  validated.  Figure  29  shows  the  percent 
error  between  the  BRLCB  computed  surface  area  and  the  Robbins  and  Lynn  surface  area.  The  maximum 
percent  error  is  approximately  0.002%. 
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Figure  25.  Deduced  burning  rate,  test  case  1. 


igure  26.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rate,  test  case  1. 


lime  (ms) 


0.022  - 
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Percent  error  between  BRLCB  generated  pressure  and  the  analytically  Renerated  pressure,  test  case  2. 


gure  29.  Percent  error  between  BRLCB  generated  sunacc  area  and  the  analytically  generated  surface  area,  test  case  3. 


The  authors  believe  that  the  results  of  these  three  test  cases  sufficiently  validate  the  computer 
implementation  of  the  mathematical  model  described  in  section  3  of  this  report  (except  for  the  ETC 
option).  Although  the  Robbins  and  Lynn  analytic  data  arc  for  a  homogeneous  single  layered  grain, 
BRLCB  utilized  the  same  algorithm  in  dealing  with  multilayered  grains  as  for  single  layer  grains.  Thus, 
the  BI^LCB  analysis  for  multilayered  grains  should  be  as  accurate  as  for  the  single  layer  gains.  The  only 
question  for  the  multilayered  grains  should  be  the  codes  ability  to  handle  the  possible  jump  discontinuities 
in  physical  properties  associated  with  going  from  one  layer  to  the  next.  At  most,  any  affect  due  to  the 
discontinuities  should  be  observed  only  in  those  time  steps  which  require  burning  in  two  distinct  layers. 
However,  as  will  be  shown  in  section  4.5,  BRLCB  is  able  to  accurately  handle  the  discontinuous  jumps 
between  grain  layers. 

4.2  Scries  II  Test  Cases.  The  three  test  cases  presented  in  Section  4.1  were  run  on  a  CRAY  computer 
with  16  significant  figures  for  the  pressure  input  data.  Unfortunately,  experimental  data  seldom  arc  tliis 
accurate.  Thus,  it  was  of  interest  to  determine  the  affect  less  accuracy  in  input  pressure  data  would  have 
on  the  computed  burning  rates.  These  results  are  presented  in  test  cases  4  and  5. 

Test  Case  4.  This  test  case  is  identical  to  test  case  1  except  that  the  accuracy  of  tlic  pressure  data  was 
reduced  from  16  significant  figures  to  6.  The  percent  error  in  the  computed  burning  rates  is  shown  in 
Figure  30. 

As  can  be  seen  in  Figure  30,  the  maximum  percent  error  is  approximately  0.03%.  This  represents  an 
order  of  magnitude  increase  in  percent  cram  over  the  0.(X)2%  observed  when  the  pressure  was  accurate 
to  16  significant  figures  (Figure  26). 

Test  Case  5.  In  this  test  case,  the  accuracy  of  the  input  pressure  data  was  further  reduced  to  four 
significant  figures.  As  shown  in  Figure  31,  the  maximum  percent  error  has  risen  to  between  3%  and  4%. 

4.3  Scries  III  Test  C.ises.  The  objective  of  the  test  cases  in  this  scries  is  to  validate  the  coding 
associated  with  grain  geometry  form  functions  (all  geometries  except  the  sandwich  grain  with  inhibited 
sides  for  which  no  pressure-time  data  could  be  genenhed).  For  test  scries  I,  the  grain  geometry  was  a 
singlc-pcrfordtcd  cylinder  and  the  results  from  that  test  scries  confirmed  that  the  form  function  was 
correctly  implemented  for  the  single-pciforaicd  grain.  Form  functions  for  11  of  the  remaining  12  grain 
geometries  available  in  BRLCB  will  be  validated  in  this  series  by  perfomimg  a  burning  rate  reduction 
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Percent  error  betv.^een  BRLCB  deduced  burning  rate  and  assumed  bumini 
rate  for  pressure  data  with  six  significant  figures  of  accuracy,  test  case  4. 


LEGEND 

100*(BR  -  BRact-ual)/BRactual 


Percent  error  between  BRLCB  deduced  burning  rate  and  assumed  burning 
rate  for  pressure  data  with  four  signiricaht  figures  of  accuracy,  test  case  !)■ 


using  BRLCB.  However,  unlike  test  series  I,  pressure-time  data  arc  not  provided  by  an  analytic  solution 
but  by  adapting  the  interior  ballistic  code  IBHVG2  (Anderson  and  Fickic  1987)  to  simulate  a 
closed-chamber  firing.  Thus,  the  pressure-time  input  data  will  contain  round-offAruncation  error  which 
may  affect  the  accuracy  of  tlie  BRLCB  bum  rate  computation  (i.e.,  increase  in  percent  error  between 
deduced  and  assumed  bum  rate).  In  addition,  tlie  IBHVG2  piessure-time  data  are  accurate  to  only  four 
decimal  places  (five  total  decimal  digits).  Therefore,  based  upon  results  of  test  scries  II  and  potential 
round-off/truncation  errors  within  the  pressure  data,  percent  cnors  between  deduced  and  assumed  bum 
rates  on  the  order  of  several  percent  should  be  expected  from  the  computations. 

Test  Case  6.  For  tlie  spherical  form  function,  the  percent  error  between  deduced  and  assumed  bum 
rate  is  shown  in  Figure  32.  Inputs  used  for  this  test  case  arc  idcnUca!  to  those  of  the  sample  test  case 
provided  in  section  2. 

Test  Cases  7-15.  Results  for  the  remaining  grain  geometries  (except  sandwich  with  inhibited  sides) 
arc  presented  in  Figures  33-42.  In  all  cases,  the  input  to  IBHVG?  and  BRLCB  remained  the  same  except 
for  the  grain  geometry.  'Hiese  fixed  inputs  are  provided  in  Table  39. 

Table  39.  Fixed  Input  Parameters  for  IBHVG2  and  BRLCB  Simulations 


Closed-chamber  volume; 

300  cm^ 

Grain  geometry: 

Varied 

Propellant: 

Impetus: 

Flame  tempctaiurc: 

Ratio  of  specific  heats: 
Molecular  weight; 
Covolume; 

Density; 

Total  mass: 

1.140  J/g 

3.410  K 

1.225 

24.869 
.996  cm^/g 

1.58  g/cnr 

90  g  +  1  g  igniter 

Bum  rate  law: 

r  =  .16  cm/s 
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Figure  32.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rate,  spherical  grain  geometry,  test  case  6. 
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Figure  33.  Percen!:  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rate,  cord  grain  eeomelrv.  test  case  7. 
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Figure  35.  Percent  error  between  BkLCB  deduced  bum  rate  and  assumed  burning  rate,  slotted  tube  grain  geometry,  test  case  9. 


i 


Figure  36.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rate,  7-perforation  grain  geometry,  tes:  case  IQ. 


Figure  37.  Perccm  error  between  BRLCB  deduced  hum  rate  and  a.ssumcd  burning  rale,  7-pcrforaiion  hex  erain  geometry,  test  case  1 1. 


Figure  38.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  raie.  19-Derforation  grain  geometry,  test  case  12. 
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Figure  39.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed 
burning  rale.  l9-pcrforation  hex  grain  geometry,  test  case  13. 
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Figure  40.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  buming  rate.  37-perforation  grain  geometry,  test  case  14. 


Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  bumiri! 
rate,  cord  with  inhibited  ends  grain  geometry,  test  case  15. 
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Figure  42.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rate,  cigarette  grain  geometry,  test  case  16. 


As  can  be  observed  from  Figures  32  through  42,  the  percent  error  from  the  BRLCB  calculations  is 
generally  within  the  2-3%  expected.  The  exceptions  being  the  19-pcrforation  hex  grain  and  the  cigarotte 
grain.  However,  as  shown  in  Figure  39  where  the  actual  data  points  are  starred,  the  percent  error  is 
outside  of  the  expected  range  for  only  two  points.  In  tlte  case  of  the  19-perforation  hex  grain,  these  points 
correspond  to  a  grain  slivering  event;  for  the  cigarette  grain,  the  points  are  the  first  two  points  in  the 
calculation  where  the  numerical  methods  utilized  will  be  most  inaccurate.  Based  upon  the  results  of  this 
test  series,  the  authors  believe  that  the  grain  geometry  fonn  functions  have  been  correctly  implemented, 

4.4  Series  IV  Test  Case.  Test  Case  17.  The  objective  of  this  series  was  to  validate  the  BRLCB  ETC 
bum  rate  reduction  option.  The  ETC  option  is  implemented  in  the  mathematical  model  by  assuming  that 
tlie  electrical  energy  is  strictly  additive  with  the  chemical  energy.  Thus,  tlie  electrical  energy  simply 
represents  a  source  term  to  the  energy  balance  equation.  As  in  the  series  III  test  cases,  input  data  were 
generated  using  an  interior  ballistic  code.  In  this  case,  a  modified  version  of  IBHVG2  (Eamhart  and 
Winsor  1992)  capable  of  simulating  the  ETC  process.  Inputs  for  the  modified  IBHVG2  simulation  and 
BRLCB  reduction  arc  provided  in  Table  40. 

Table  40.  Input  Para.rietcrs  for  IBHVG2  and  BRLCB  Simulations  for  ETC  Bum  Rate  Reduction 


Closed-chamber  volume: 

65.5  cm^ 

Grain  geometry: 

Sphere 

Propellant: 

Impetus: 

1,165  J/g 

Flame  temperature: 

3,680  K 

Ratio  of  specific  heats: 

1.2155 

Molecular  weight: 

26.2532 

Covoluinc: 

.996  cra^/g 

Density: 

1 .6  g/cm^ 

Total  mass: 

13  g  +  .5  g  igniter 

Bum  rate  law: 

r  =  .2259  cm/s 

The  electrical  energy  added  to  the  system  is  shown  in  Figure  43. 


Percent  error  between  deduced  and  assumed  bum  rate  is  given  in  Figure  44.  As  shown  in  the  figure, 
the  percent  error  is  well  within  the  expected  limits.  Thus,  the  authors  feel  the  ETC  bum  rate  reduction 
option  is  properly  implemented. 

4,5  Series  V  Test  Cases.  In  this  series,  the  ability  of  BRLCB  to  handle  layered  propellants  with 
constant  and  varying  properties  is  investigated.  Test  case  18  addresses  constant  properties,  while  varying 
properties  are  investigated  in  test  case  19.  Unlike  the  previous  test  cases,  these  two  test  cases  use  the 
pressure  generation  option  of  BRLCB  to  generate  the  pressure  data  used  in  the  bum  rate  reduction. 
Although  this  may  appear  to  represent  a  circular  validation,  this  is  not  the  case.  The  object  of  these  test 
cases  is  to  determine  if  BRLCB  can  cleanly  handle  the  potential  discontinuous  jumps  in  properties 
(thermochemical  and  bum  rate)  at  the  layer  boundaries.  The  pressure  generation  option  of  BRLCB 
produces  a  pressure  history  at  equal  time  steps,  the  jump  between  layers  docs  not  necessarily  occur  at  one 
of  these  time  steps.  Therefore,  the  pressure  history  generated  by  BRLCB  "smears"  the  jump  between 
layers  across  one  time  step  and  provides  no  additional  information  to  the  bum  rate  reduction  option  that 
would  not  be  available  from  an  experimental  closed-chamber  firing. 

Test  Case  18.  In  this  test  case,  a  five-layer  sphere  with  constant  properties  in  each  layer  is 
investigated.  Input  information  is  provided  in  Table  41. 

The  assumed  bum  rates  used  in  this  case  arc  shown  in  Figure  45.  The  corresponding  pressure-time 
data  generated  by  BRLCB  are  shown  in  Figure  46. 

Results  of  the  BRLCB  calculation  arc  presented  in  Figures  47  and  48.  Figure  47  is  the  deduced  bum 
rate  and  Figure  48  the  percent  error  between  the  deduced  and  calculated  burning  rates. 

As  illustrated  in  Figure  47,  BRLCB  is  able  to  cleanly  handle  the  jumps  between  the  layers.  The 
percent  error  shown  in  Figure  48  is  felt  acceptable  by  the  autiiors,  being  less  than  0.01%  for  all  but  the 
last  10  points 

Test  Ca.se  19.  This  test  case  is  similar  to  test  case  18  except  that  the  properties  in  each  layer  (except 
bum  rate)  arc  varying  as  shown  in  Table  42. 


78 


Electrical  energy  added  for  ETC  ciosed-chamber  1BHVG2  and  BRLCB  simulations,  test  case  17. 


Figure  44.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rate,  test  case  17. 


Figure  45.  Assumed  bum  rates  for  five-layer  sphere  with  constant  properties  in  each  layer,  test  case  18. 


igure  47.  Deduced  burning  rates  for  flve-Iayer  sphere  with  constant  properties  in  each  layer,  test  case  18. 
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Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rale 
five-layer  sphere  with  constant  properties  in  each  layer,  test  case  1 8. 


Tabic  41.  Input  Parameters  for  BRLCB  Simulations  for  Five-Layer  Sphere  With 
Constant  Properties  in  Each  Layer,  Test  Case  1 8 


Closed-chamber  volume:  3CX)  cm^ 


Grain  geometiy:  sphere  (five  layers)  (.5  cm  diameter) 


Starting  depth: 

(cm) 

Layer  1 

0.0 

Layer  2 

.0263 

Layer  3 

.06688 

Layer  4 

.1316 

Layer  5 

.23434 

Propellant  properties: 

Impetus  (J/g): 

382.5 

540.4 

803.7 

1,064.2 

1,398.3 

Flame  temperature  (K): 

2,300 

2,600 

2,900 

3,200 

3,700 

Gamma  (-): 

1.27 

1.26 

1.25 

1.24 

1.22 

Molecular  weight  (-): 

50 

40 

30 

25 

22 

Covolume  (cm^/g; 

.95 

.92 

.89 

.85 

.80 

Density  (g/cm^): 

1.4 

1.45 

1.5 

1.55 

1.65 

Total  mass:  1 10  g  -t-  2  g  igniter 

Bum  rate  law 

(cm/s): 

r=1.0P' 

r=.02P'-^ 

rt=.2P^ 

r+.OlP''^ 

r=.05P'-'’ 

The  assumed  bum  rates  used  in  this  case  arc  shown  m  Figure  49.  The  corresponding  prcssurc-time 
data  generated  by  BRLCB  arc  shown  in  Figure  50  by  the  solid  line. 

Results  of  the  BRLCB  calculation  arc  prc.scntcd  in  Figures  51  and  52.  Figure  51  is  the  deduced  bum 
rate  and  Figure  52  the  percent  error  between  the  deduced  and  calculated  burning  rates. 

As  with  the  test  case  18,  Figure  51  indicates  that  BRLCB  is  capable  of  cleanly  handling  the 
discontinuous  jumps  between  layers.  The  percent  error  shown  in  Figure  52  is  almost  identical  to  that 
shown  in  Figure  48  for  the  5  layer  sphere  with  constant  properties.  Thus,  the  use  of  the  integrals  to 
compute  tire  average  of  the  variable  properties  appears  to  have  litUc  or  no  effect  on  tlic  accuracy  of  the 
bum  rate  calculation  compared  to  having  constant  property  values. 
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Table  42.  Input  Parameters  for  BRLCB  Simulations  for  Five-Layer  Sphere 
With  Varj'ins  Properties  in  Each  Layer,  Test  Case  19 


Qosed-chambcr  volume:  300  cm^ 

Grain  geometry:  sphere  (five  layers)  (.5  cm  diameter) 

Starting  depth: 

Layer  1 

Layer  2 

Layer  3 

Layer  4 

Layer  5 

(cm) 

0.0 

.02636 

.06688 

.1316 

.23434 

Propellant  piopcnies: 

Impetus  (J/g): 

Start: 

382.5 

540.4 

803.7 

1,064.2 

1,398.3 

End: 

Flame  temperature  (K): 

540.4 

803.7 

1,064.2 

1,398.3 

1,398.3 

Start: 

2,300 

2,500 

2,900 

3,200 

3,700 

End: 

Gamma  (-): 

2,600 

2,900 

3,200 

3,700 

3700 

Start: 

1.27 

1.26 

1.25 

1.24 

1.22 

End: 

Molecular  weight  (-): 

1.26 

1.25 

1.24 

1,22 

1.22 

Start: 

50 

40 

30 

25 

22 

End: 

40 

30 

25 

22 

22 

Covolumc  (cm^/g): 

Start: 

.95 

.92 

.89 

.85 

.80 

End: 

.92 

.89 

.85 

80 

.80 

Density  (g/cm^): 

Start: 

1.4 

1.45 

1.5 

1.55 

1.65 

End: 

1.45 

1.5 

1.55 

1.65 

1.65 

1  Total  mass:  110  g  -i-  2  g  igniter  | 

Bum  rate  law 

(cm/s): 

4  =  1.0p5 

r  =  .02P*-^ 

r  =  .2P’ 

r  =  .OIP^-' 

r  =  .05P^-^ 
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igurs  52.  Percent  error  between  BRLCB  deduced  bum  rate  and  assumed  burning  rate, 
five-layer  sphere  with  varying  properties  in  each  layer,  test  case  19. 


5.  COMPARISON  WITH  OTHER  CLOSED-CHAMBER  DATA  REDUCTION  PROGRAMS 


In  this  section,  bum  rate  reduction  results  from  BRLCB  will  be  discussed  in  comparison  to  results 
from  other  closed-chamber  data  reduction  codes.  The  objective  is  not  validation  of  the  mathematical 
model  and  computer  implementation,  but  rather  to  provide  the  user  with  a  reference  point  for  use  when 
comparing  or  utilizing  results  from  different  reduction  programs.  In  addition  to  BRLCB,  three  other  data 
reduction  programs  have  been  or  are  still  being  employed  by  the  Weapons  Technology  Directorate,  ARL, 
to  reduce  closed-chamber  data.  These  programs  are  CBRED2  (Juhasz  and  Price  1977);  MINICB  (Oberle, 
Juhasz,  and  Griffie  1987);  and  SIMPCB  (Robbins  and  Horst  1976).  In  a  comparison  of  the  results 
produced  by  SIMPCB  and  BRLCB  for  JA2  stick  propellant,  Robbins  (1991)  found  virtually  no  difference 
in  the  computed  bum  rate.  Thus,  comparison  with  SIMPCB  will  not  be  included  in  tliis  rep>ort.  In  an 
earlier  report  (Oberle,  Juhasz,  and  Griffie  1987),  it  w.as  shown  that  CBRED2  and  MINICB  calculate  bum 
rate  equal  to  about  1%.  Since  CBRED2  is  no  longer  available  at  the  Weapons  Technology  Directorate, 
ARL,  results  from  BRLCB  will  be  compared  to  MI^^CB  with  the  relation  to  CBRED2  inferred  from  die 
earlier  comparison  of  CBRED2  and  MINICB. 

For  the  comparison,  a  19-perforation  grain  was  selected  since  this  was  the  same  grain  geometry  used 
in  the  comparison  of  MINICB  and  CBRED2.  Details  of  the  closed-chamber  experiment  arc  provided  in 
Table  43. 


Table  43.  Qosed  Chamber  Parameters 


.  .  ■  ■■  -  - —  , 

Igniter:  1.5-g  black  powder 

Propellant;  72.2715  g 

Geomeuy:  10-pcrforation 

Thermochemistry 

1 

Length  (cm) 

Diameter  (cm) 

Perforation  diameter  (cm) 

Inner,  outer,  and  middle  web  (cm) 

0.722884 

0.865886 

0.035560 

0.114681 

Impetus  (J/g) 

Flame  temperature  (K) 
Density  (g/cm^) 
Molecular  weight  (-) 
Covolumc  (cmr/g) 
Gamma  (-) 

1,069,0 

2,671.0 

1.636 

20.782 

1.165 

1.2688 

Bomb:  211  cm^ 

A  comparison  of  the  computed  bum  rate  is  provided  in  Figures  53  and  54.  Figure  53  provides  the 
complete,  unfiltered,  nondccimated  output  from  BRLCB  and  the  decimated  (data  points  deleted)  output 
from  MINICB  (MINICB  automatically  decimates  the  output  to  reduce  the  data  to  be  processed). 
Figure  54  is  filtered  (low  pass  filter)  BRLCB  bum  rate  data  vs.  the  MINICB  results.  A  comparison  of 
the  bum  rate  in  the  figures  indicates  very  little  difference  in  the  computed  bum  rate  except  below  10  MPa 
and  above  approximately  300  MPa.  The  differences  above  300  MPa  can  be  attributed  to  the  difference 
in  the  form  function  used  to  describe  the  grain  geometry  after  "silvering."  BRLCB  explicitly  accounts 
for  "silvering"  whereas  MINICB  does  not.  This  results  in  the  computed  bum  rate  for  BRLCB  remaining 
near  linear  in  the  300-400  MPa  range  (Figure  54).  To  quantify  the  results.  Table  44  compares  the  bum 
rate  results  for  selected  pressure. 


As  illustrated  in  Figure  53  and  54  and  Table  44,  BRLCB  and  MINICB  produce  bum  rate  rcsulLs  to 
witliin  1%  ( 100-350  MPa)  and  3%  at  50  MPa.  Considering  that  differences  in  the  computational  schemes 
and  grain  geometiy  foim  functions,  the  authors  consider  the  results  to  be  equivalent.  Tlius,  BRLCB,  at 
least  for  the  limited  number  of  cases  investigated,  produces  equivalent  results  when  compared  to  other 
closed-chamber  data  reduction  programs  used  at  the  Weapons  Technology  Directorate,  ARL. 


Table  44.  Comparison  of  Bum  Rate  BRLCB  vs.  MINICB 


Pressure 

(MPa) 

Bum  Rate  -  BRLCB 
(cm/s) 

Bum  Rate  -  MINICB 
(cm/s) 

%  Difference 
(-) 

50 

3.060 

3.152 

-t-3.00 

100 

5.335 

5.390 

-1.03 

150 

7.168 

7.165 

-0.04 

200 

9.941 

10.000 

0.60 

250 

13.442 

13.340 

0.73 

300 

17.325 

17.280 

-0.26 

350 

20.300 

20.220 

-0.39 

400 

21.217 

20.900 

-1.50 

6.  CONCLUSIONS 


BRLCB  is  a  closed-chamber  dat  i  analysis  progrturi  designed  to  accommodate  homogeneous,  layered, 
and  deterred  propellants.  ITic  program  offers  five  modes  of  operation:  1)  bum  rate  analysis,  2)  synthetic 
pressure-time  history  generation,  3)  surface  area  analysis,  4)  intcmipicd  burner  analy.sis  and  5)  ETf'C  bum 
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rate  reduction.  The  incorporation  of  continuously  varying  thermochcniical  properties  allows  aU  propellant 
geometries  to  be  treated  as  layered  or  deterred.  Comparisons  with  other  closed-chamber  data  analysis 
programs  indicates  essentially  identical  results. 


-OL 


larison  of  bum  rate  results;  BRLCB  vs.  MINICB 
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APPENDIX  A: 


REQUIRED  INFORMATION  FOR  RUNNING  BRLCB 


INTENTIONALLY  LEFT  BLANK. 
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Option  1:  Create  Master  Information  File 

A.  General  -  name  of  master  information  file  if  current  file  being  created  from  previous  file 

B.  Identification  information 

1 .  Project  name 

2.  Person  requesting  work 

3.  Master  information  file  name  (.MAS  extension  suggested) 

C.  Propellant  information 

1 .  Propellant  type 

2.  Propellant  lot 

3.  Propellant  source 

4.  Number  of  layers  in  grain 

5.  Constant  or  varying  grain  propcnics 

6.  Flame  temperature  (K)  for  each  layer 

7.  Impetus  (J/g)  for  each  layer 

8.  Molecular  weight  (-)  for  each  layer 

9.  Covolume  (cc/g)  for  each  layer 

10.  Ratio  of  specific  heats  (-)  for  each  layer 

11.  Propellant  density  (g/ec)  for  each  layer 

12.  Propellant  grain  geometry 

13.  Grain  diameter  (cm) 

14.  Grain  length  (cm)  if  needed 

15.  Perf  diameter  (cm)  if  needed 

16.  Inner  web  (cm)  if  needed 

17.  Middle  web  (cm)  if  needed 

18.  Outer  web  (cm)  if  needed 

D.  Igniter  information 

1.  Igniter  name 

2.  Igniter  lot 

3.  Ingitcr  source 

4.  Igniter  impetus  (J/g) 

5.  Igniter  flame  temperature  (K) 

6.  Igniter  density  (g/cc) 

7.  Igniter  molecular  wciglit  (-) 

8.  Igniter  covolume  (cc/g) 

9.  Igniter  ratio  of  specific  heats  (-) 

Option  2:  Update  gage  information 

1.  Gage  ID 

2.  Gage  calibration  date 

3.  Gage  calibration  coefficients  (.sc  ond  order  fit  to  produce  MPa) 
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Option  3:  Prepare  pressure/time  data 

1.  Foimat  of  pressure/time  data  (see  Table  13) 

2.  Name  of  pressure/time  data  file 

3.  Output  file  name  for  pressure/time  data  (all  files  for  the  reduction  will  use  this  file  name) 

4.  Time  step  for  data  in  milliseconds 

5.  Using  VuPoint/BRL  procedure  -  calibration  voltage  of  charge  amplifier 
Option  4;  Repare  firing  information  file 

1.  Name  of  master  file  to  be  utilized  in  creating  the  firing  information  file 

2.  Type  of  computation  (sec  table  17) 

3.  Name  of  pressure/time  file  created  in  Option  3  to  be  used  in  tlie  reduction 

4.  Closed  chamber  volume  (cc) 

5.  Initial  temperature  of  closed  chamber  (K) 

6.  Propellant  mass  (g) 

7.  Ignition  mass  (g) 

8.  Initial  propellant  temperature  (K) 

9.  Initial  igniter  temperature  (K) 

10.  Ranges  for  bum  rate  low  computations 

11.  For  El'C  option  -  name  of  electrical  energy  file 

Option  5:  Smooth  pressure/time  data 

1.  Type  primer 

2.  File  name  (same  name  as  firing  infoimation  file) 

3.  Bridge  lengtli  option 

4.  Number  of  wildpoint  removal  pa.s.scs 

5.  Wildpoint  tolerance 

6.  Number  of  smoothing  passes 

7.  Bridge  length  (fixed  bridge  length  option  only) 

8.  Pressure  in  Ml^a  associated  with  grain  slivering  (if  slivering  option  selected) 

9.  Make  pressure  strictly  monotoiiic?  (yes/no) 

10.  Vivacity  information  required?  (yes/no) 

Option  6;  Perform  data  analysis 

1.  Set  closed  ehainber  air  mass  to  0.0?  (yes/no) 

2.  Change  heat  loss  fraction?  (yes/no)  If  yes,  new  value  for  heat  loss  fraction. 

Option  7;  Prepare  output 

1.  File  name 

2.  Date  of  closed  chamber  firing 

3.  Comments  on  firing 

4.  Comments  on  reduction 

5.  Type  output  (hard  copy,  file  or  both) 

6.  Printer  port  (LPTl,  etc.) 

7.  Number  of  lines  per  primed  page 
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8.  Skip  factor  for  printing  tabular  data 

9.  FFT  on  bum  rats?  (yesAio) 

10.  Printer  type 

11.  Pressure  ranges  for  bP^  bum  rate  laws 

12.  Summary  output  sheet?  (yesAio) 
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INTENTIONALLY  LEFT  BLANK. 
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APPENDIX  B: 

FILE  STOUCTURE:  .MAS  AND  .INF  PILES 


105 


■MSB 


Lines  1-6:  Array  A2(l)  through  A2(6) 

Format  A80 

Comments  on  the  firing  and  reduction 

Lines  7-26:  Array  A  1(1)  through  A  1(20) 
Format  A20 


Line  7:  A  1(1) 

Project  name 

Line  8:  Al(2) 

Person  requesting  work 

Line  9:  A  1(3) 

Master  hie  name 

Line  10:  Al(4) 

INF  flic  name 

Line  11:  Al(5) 

.PVT  file  name 

Line  12:  Al(6) 

Propellant  type 

Line  13:  Al(7) 

Propellant  lot 

Line  14:  Al(8) 

Propellant  source 

Line  15:  A  1(9) 

Name  of  grain  geometry 

Line  16:  Al(lO) 

Output  hie  name 

Line  17:  Al(ll) 

Igniter  name 

Line  18:  Al(12) 

Igniter  lot 

Line  19:  Al(13) 

Igniter  source 

Line  20:  Al(14) 

.PDT  hie  name 

Line  21:  Al(15) 

Gage  identiheation  number 

Line  22:  Al(16) 

Bomb  type  (closed  or  interrupted) 

Line  23:  Al(17) 

Graphics  hie  name 

Line  24:  Al(18) 

Computed  bum  rate  data 

Line  25:  Al(19) 

Electrical  energy  hie  name 

Line  26:  Al(20) 

Date  of  firing 

Lines  27  through  126:  Array  A3(l>-A3(100) 


Line  27:A3(l) 
Line  28:A3(2) 

Line  29:  A3(3) 
Line  30:  A3(4) 
Line  31:  A3(5) 

Line  32:  A3(6) 
Line  33:  A3(7) 
Line  34:  A3(8) 
Line  35:  A3(9) 
Line  36:  A3(10) 
Line  37:  A3(ll) 
Line  38:  A3(12) 
Line  39:  A3(13) 
Line  40:  A3(14) 
Line  41:  A3(15) 
Line  42:  A3(16) 
Line  43:  A3(17) 


Maximum  grain  depth  for  burning 

Calculation  mode:  (l=bum  rate;  2=prcssure  generation;  3-surfacc  area 
4=intcmiptcd  burner.  5=ETC  bum  rate) 

Mass  (grams)  of  a  single  grain 
Number  of  layers  per  grain 

Number  of  bum  rate  pairs  (negative  number:  bP^n  laws; 
positive  number:  tabular  pressure-bum  rate  information) 

Not  used 

Grain  length  (cm) 

Grain  outer  diameter  (cm) 

Grain  perf  diameter  (cm) 

Grain  inner  web  (cm) 

Grain  middle  web  (cm) 

Grain  outer  web  (cm) 

Igniter  impetus  (J/g) 

Igniter  flame  temperature  (K) 

Igniter  density  (g/cc) 

Igniter  molecular  weight 
Igniter  covolumc  (cc/g) 
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Line  44:  A3(18)  Igniter  ratio  of  specific  heats 

Line  45:  A3(19)  Gage  input  voltage  (volts) 

Line  46:  A3(20)  Constant  in  gage  fit 

Line  47:  A3(21)  Linear  term  in  gage  fit 

Line  48:  A3(22)  Quadratic  in  gage  fit 

Line  49:  A3(23)  Bomb  volume  (cc) 

Line  50:  A3(24)  Initial  bomb  temperature  (K) 

Line  51:  A3(25)  Total  propellant  mass  (grams) 

Line  52:  A3(26)  Igniter  mass  (grams) 

Line  53:  A3(27)  Initial  propellant  temperature  (K) 

Line  54:  A3(28)  Initial  igniter  temperature  (K) 

Line  55:  A3(29)  Pressure  due  to  igniter  (MPa) 

Line  56:  A3(30)  Total  electrical  energy  (MJ) 

Line  57:  A3(31)  Flag  for  grain  type  (()=variable  layers:  l=homogencous,  one  layer, 

2=homogeneous  layers,  several  layers) 

Line  58:  A3(32)  through  Line  67:  A3(41)  Not  used 
Line  68:  A3(42)  Time  step  in  milliseconds 
Line  69:  A3(43)  Observed  maximum  pressure 
Line  70:  A3(44)  Theoretical  maximum  pressure 
Line  71:  A3(45)  Grain  geometry  number 
Line  72:  A3(46)  Number  of  propellant  grains 

Line  73:  A3(47)  Not  used 

Line  74:  A3(48)  Not  used 

Line  75:  A3(49)  Not  used 

Line  76:  A3(50)  tltrough  Line  125:  A3(99)  contain  the  bum  rate  law  information;  absolute  value  of  A3(5) 
determines  the  number  of  pairs.  If  A3(5)  is  negative,  then  bP^n  bum  rate  laws  arc 
used.  These  laws  are  stored  in  the  following  manner. 

A3(50)  -  first  law  coefficient,  A3(51)  -  corresponding  exponent, 

A3(52)  -  second  law  coefficient,  A3(53)  -  corresponding  exponents,  etc. 

If  A3(5)  is  positive,  the  tabular  pressure-bum  rate  infonnalion  is  stored  in  the 
following  manner. 

A3(50)-A3(79):  Pressures  in  MPA 
A3(80)-A3(99):  Corresponding  bum  rate  (cm/s) 

Line  126:  A3(100)  Number  of  computational  steps  performed 

The  remaining  825  lines  of  the  file  contain  properties  of  a  single  grain  of  the  propellant  and  storage 
locations  for  the  computation.  These  825  lines  arc  thought  of  as  a  three-dimensional  array,  PX(1 1,15,5). 
The  information  is  written  in  the  following  manner. 

DO  10  1=1,11 
DO  20  J=l,15 
DO  30  K=l,5 
WRITE(-,*)PX(U,K) 
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30 

CONTINUE 

20 

CONimJE 

10 

CONTINUE 

The  stnicturc  of  the  array  PX  is  as  follows: 

PXa.J.K): 

I=pnoperty 

1  =  Depths  into  grain  (0.0  is  depth  of  outer  layer) 

2  =  Impetus  (J/g) 

3  =  Flame  temperature  (K) 

4  =  Density  (g/cc) 

5  =  Molecular  weight  (-) 

6  =  Covolume  (cc/g) 

7  =  Ratio  of  specific  heats,  gamma 

8  =  Volume  if  K  is  3,  4,  or  5 

Value  of  1  for  K  a  i  or  2 

9  =  Universal  gas  constant  (cal/g-K) 

10  =  Ratio  cf  specific  heat,  constant  volume  (cal/g-K) 

1 1  =  Specific  energy  (cal/g) 

J  =  layer 

K  =  determines  specific  content  of  array  cell 

1  =  value  at  beginning  of  layer 

2  =  value  at  end  of  layer 

3  =  value  of  integration  over  entire  layer 

Integral  =  Integral  (Property  *  Density  *  dx) 

4  =  value  of  last  converged  integral 

5  =  value  of  last  computed  integral 

Exceptions: 

PX(1,-,1):  beginning  depth  of  layer 

PX(l,-,2):  depth  or  len^  of  layer,  not  the  depth  at  the  end  of  the  layer 

PX(l,-,3):  Not  used 

PX(l,-,4):  depth  for  la.st  converged  integral 

PX(l,-,5):  depth  for  last  computed  integral 
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iNTENnONALLY  LEFT  BLANK. 
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APPENDIX  C: 
SAMPLE  OUTPUT 


Ill 


INIT-NTIONALLY  LEFT  BLANK. 


BURNING  RATE  ANALYSIS 
BRLCB  V3.0 

ADVANCED  BALLISTIC  CONCEPTS  BRANCH  -  BRL, 


Project;  Test  of  7-perf 
Inf  File:  a;pcrf7.inf 
P/T  File:  a:perf7.pvt 
Smoothed:  a:perf7.pdt 
Fired  on:  2/23/92 


Requested  by:  obcrle 

Created  From  .MAS  File:  a;pcrf7.mas 

Calculation  Output  File:  a:perf7.out 

Graphics  File:  a:pcrf7.dat 


IGNFfER  INFORMATION 

The  igniter  used  is:  Dummy  JA2  Lot:  Dummy 

The  source  for  the  igniter  is:  William  Obcrle 

IGNITER  THERMOCHEMICAL  PROPERTIES 
Impetus  (J/g):  1,140.0  Molecular  Weight:  24.86900 

Flame  Temperature  (K):  3,410.0  Covolume  (cc/g):  .99600 

Density  (g/cc):  1.58000  Gamma:  1.22500 

PROPELLANT  INFORMATION 
The  propellant  used  is:  JA2  Lot;  Dummy 

The  source  for  tlic  propellant  is:  obcrle 

Propellant  Thcnmochcmical  Properties:  Following  Sheets  of  Output 

PROPELLANT  GRAIN  GEOMETRY 

Grain  Type:  7-Pcr.  Cyl. 

Length  (cm):  2,0000(X) 

Outer  Diam.  (cm):  1.3(X)000 

Perf  Diam.  (cm):  .100000 

Inner  Web  (cm):  .200000 

Outer  Web  (cm);  .30(X)00 

Bomb  Information _  Gage  Information 

Bomb  Type:  Closed  Chantber  Gage  ID:  None 

Bomb  Vol  (cc):  300.0  Input  Voltage:  .0000 

Constants  for  Fit:  A4-Bx-»-C'^2 
A:  .OOOOOE-K)0 
B:  .OOOOOE+00 
C:  .OOOOOE+00 


Temtyraturc  and  Charge  Ma.ss  Information 


Propellant  Mass  (g):  90.0000 

Initial  Temp.  Prop.  (K):  294. 

Initial  Bomb  Temp.  (K):  294. 

Number  of  Propellant  grains;  22.38 


Igniter  Mass  (g);  1 .0000 

Igniter  Temp.  (K):  294. 
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Total  #  Layers  =  1 

Chamber  Volume  (cm3)  =  300.000 
Hcat-Loss-Fraction  (n-d)  =  .000 

Time  Step  (ms)  =  .500000(X)E-01  Max  Time  steps  =  1,200 

Convergence  Criterion  =  .lOOOOOOOE-04 

Computation  type  (]=BURN-RATE.  2~PTGEN.  3=SURF  AREA,  4=INrERR  BOMB,  5=E'rC) 

1 

*  *  Ik  *  4i  111  4<  *  «  *  4i  *  41 4t  %  4<  %  *  4<  >)<  41  *  *  «  4<  4<  *  41  *  41  Ik  4(>|i  i(>  4>  %  A  *  4i  <•  4i  <ti  4>  *  %  %  *  i|<  ifr  <l>  i«i  f  *  *  i4i  <ti  If  %  %  >t>  <•<  4<  *  *  %  >l<  4<  >l<  i|< 

Bcginning(l)  of  Layer  End(2)  of  Layer 

f  f  f  4i  4i  4i  4i  f  4>  4i  4i  %  f  If  4i  4<  If  4i  f  If  If  %  If  If  4.  if  4i  4i  4<  4i  f  f  If  4>  f  f  4i  4i  4i  4i  f  4i  4i  4i  4i  *  4i  4i  4: 4i  4i  f  If  4i  4i  4t  4i  f  4i  4i  4i  4i  f  If  If  If  4i  4i  f  4i  If  If  4i  If  If  *  4>  4i  4i  f  If  f  4i  4<  4i 

Propellant  Mass  in  grams 

4i  f  4i  4i  4i  4>  4<  4i  4i  4i  4i  4i  4i  4>  4t  4i  4i  4i  f  4i  4i  4i  4i  f  4i  If  4i  4i  4i  4i  4i  4i  4i  4i  4i  4i  4i  4i  4i  4i  4>  4i  *  *  4i  4i  If  f  *  v  4<  f  4<  4i  4<  4i  4i  f  If  4<  4i  If  4i  4<  If  4<  4i  4i  If  f  4i  If  4i  4>  If  f  4<  If  If  4>  f  If  4i  4i  4i 

1  90.00002 

Igniter  1.00000 

Air  .00000 

4i  f  4i  4i  f  4i  4i  4i  f  4i  4i  4i  f  4i  If  4i  4<  4i  4i  4i  4i  4i  4i  4i  4i  If  If  4i  4i  f  f  f  f  4<  If  4<  f  If  4<  4i  f  4i  If  4i  4i  4i  4i  If  4>  4i  4i  4i  4i  4i  4i  4i  4i  If  f  4i  4i  4i  4i  4>  4i  If  4i  If  4i  4i  If  4i  4i  4i  If  f  4i  4c  4i  If  4i  f  4. 4i  4i 

Propellant  Density  in  g/cc 

4i  4i  4i  If  4  4>  4i  4i  4i  4i  4i  4i  4i  4<  4i  4i  If  4i  f  4i  4i  4i  4i  4i  4i  4: 4i  4i  4i  4>  4<  4i  4c  4<  4i  41  If  f  4i  4i  4i  41 4i  4i  4i4i  4i  f  4i  4i  4>  If  4>  If  If  4>  w  f  If  If  4<  4i  4i  f  If  If  4i  4c  4i  4i  4<  4i  4<  4i  4i  4c  If  4<  If  4i  f  4' 4c  4i  4i 

1  1.58000  1.58000 

Igniter  1.58000 

4i  4i  4i  4i  4i  4i  4i  4i  4i  4i  4c  4c  4i  4i  4i  4i  4c  4c  4c  4i  4c  4c  4i  41 4i  f  4i  4c  4i  4i  4i  4c  4c  4i  4i  4i  f  4c  4c  If  4i  4i  4i  41 4i  4i  4i  4i  If  4c  4c  4i  4i  4i  4  4i  4i  4i  4i  4i  If  4i  4i  4i  4i  4i  4i  4c  4i  4i  4i  4i  4>  4i  4i  4i  4i  4i  4c  If  If  4i  4i  4i  4i 

Propellant  Impetus  in  J/g 

444444444if4444444444444444444444444444444444:f4444444444444444444if  44444444444444444444 

1  1140.00000  1140.00000 

Igniter  1 140.00000 

4444444444444444444444444444444444444444444444444444444444444444444444444444444444444 

Propellant  Flame  Temperature  in  deg  K 

4444444444444444444*44444444444444444444444444444444444444444444444444444444444444444 

1  3410.00000  3410.00000 

Igniter  3410.00000 

444444444444444444444444.444444444444*444444444444444444444444444444444444444444*44444 

Propellant  Molecular  Weight  in  g/g-molc 

4444444444444444444444444444444444444444444444444444444444444444444444444444444*44444 

1  24.86907  24.86907 

Igniter  24.86900 

Air  28.90000 

44444444444444444444444444444444444444*4444444444444444444444444444444444444444*44444 

Propellant  Gas  Constant  in  cal/g-deg  K 

44444444444444444*44444444444444444444444444444444444444*444444444444444444444444*444 

1  .07991  .07991 

Igniter  .07991 

Air  .06876 

44444444444444444444444444444444444444444444444*4444444444444444444444444444444444444 

Gamma  (Ratio  of  Specific  HcaLs) 

4c  4i  4(  94(  %  4c  %  4(  4(  9^  4t  4i  %  4c 4i  4;  3t(  ih  %  ift  4<  %  4(  4*  %  ^  4(  %  %  %  %  %  >4(  >4i  4(  4(  %  4(  4( 

1  1.22500  1.22500 

Igniter  1.22500 

Air  1.40000 

444444444444444444444444444444444444444444444*44444444444444444444444444444444444*444 
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Specific  Heat  (Cv)  in  cal/g-dcg  k 

41 4i  4(  4i  4(  *  *  *  4i  *  4  *  4(  4i  1(1 4 4r  %  4<  *  4c  4i  lit  4<  4<  «  4  If  itr  !<(  «  4i  )|i  4r  Ik  *  4(  #1 >*1 *  4i  4i  4(  %  4i  4i  %  4i  4> 

1  .35513  ,35513 

Igniter  .35514 

Air  .17190 

4i  4i  4i  4i  4i  4(  4(  4i  4i  4i  if  4i  4<  4c  4c  4>  4c  4(  4c  4i  ill  >ii  4<  4t  4c  4c  4c  4c  4(  4i  4t  4i  4c  4(  4i  4i  4(  4i  4c  4<  4(  4c  4i  4t  4c  4t  4(  4t  4(  4(  4(  4(  4c  4c  4t  4|  4i  4i  4(  4(  4|  4i  4i  4i  4i  4i  4i  4i  4i  4c  4i  4|  4|  4|  4>  4>  4c  4c  4t  4<  4|  4t  4>  4c 

Covolume  in  cm3/g 

4t  4i  4<  4c  4(  4c  4c  4c  4t  4c  4c  4c  4c  4>  4c  4c  4t  4c  4c  If  4>  4c  4c  4<  4c  4<  4t  4c  4>  4c  4c  4c  4c  4c  4c  4i  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4t  4c  4c  4c  4>  4c  4c  4c  4c  4c  4i  4(  4t  4i  4i  4c  4c  4c  4c  4i  4(  4c  4i  4c  4c  4t  4c  4c  4c  4c  4c  4c  4c 

1  .99600  .99600 

Igniter  .99600 

Air  .98000 

4i 4i 4c  4i 4c 4c  4t 4c 4c 4c  4c 4c 4c 4c 4i  4c  4c 4c 4c 4c  4c 4c 4c  4t4i 4c 4c  4c 4c 4cif 4c 4c 4c 4c 4  4c  4c 4  f  4  4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4:4c4c4c4c4c4c4c4c4c4i4c4c4c4c4c4c4c4c4c4c4c4‘4c4c4c4c4c4c 
4c  4c  4<  4c  4c  4c  4c  4c  4c  4c  4c  4c  4t  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4i  4i  4c  4i  4c  4c  4c  4c  4i  4c  If  4  4c  4c  4c  4c  4t  *  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c  4t  4c  4c  4c  4c  4c  4i  4c  If  4c  4c  4c  4c  4c  4i  4t  4c  4c  4i  4<  4c  4c  4i  4c  4c  4c  4c  4c  4c  4c  4c  4c  4c 


******+  Maximum  Chamber  ******* 

Theoretical  Maximum  Chamber  Pressure  (MPa)  -• 
Theoretical  Maximum  Chamber  Temperature  (deg  K)  = 
***  Assuming  a  Hcat-Loss-Fraction  = 

Maximum  Total  Heat  Loss  From  Qiambcr  (cal)  = 
Maximum  Chamber  Pressure  (MPa)  w/Hcat  Loss  = 
Maximum  Qiambcr  Temperature  (deg  K.)  w/Heat  Loss  = 


.49552060E+03 

.34100000E+04 

.000 

-.43689020E-02 

.49552060E+03 

.34100000E+04 


**  Results:  Deduced  Burning  Rates  Based  on  a  Given  Chamber  Pressure  vs.  Time  **' 

***  N  ***  LAYER  ***  TIME(ms)  ***  PCH(MPa)  ***  TPMR(gm)  ***  RBR(cm/s)  “*  ASUR(cm2)  *** 


TCH(dcgK)  ***  DEPTH(cm)  *■“* 
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Pressure 

Rate 

Pressure 

8.000 

1.048 

172.000 

12.000 

1.529 

176.000 

16.000 

1.982 

180.000 

20.000 

2.450 

184.000 

4.000 

2.909 

188.000 

28.000 

3.336 

192.000 

32.000 

3.776 

196.000 

36.000 

4.206 
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40.000 

4.559 

204.000 
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208.000 
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212.000 
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5.805 

216.000 
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232.000 
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240.000 
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248.000 
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252.000 
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256.000 
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272.000 
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11.997 

280,000 
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12.364 

284.000 
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12.727 

288.000 
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13.126 

292.000 

132.000 

13.472 
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136.000 

13.844 
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148.000 

14.972 
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152.000 

15.312 
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156.000 

15.686 
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160.000 

16.052 

324.000 

164.000 

16.424 

328.000 

168.000 

16.779 

332.000 

Rate 

Pressure 

Rate 

17.158 
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31.514 
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340.000 

31.868 

17.871 

344.000 

32.194 

18.241 
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18.584 

352.000 

32.825 

18.931 

356.000 

33.173 

19.310 
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33.525 

19.658 
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33.863 

20.022 

368.000 

34.216 

20.372 

372.000 

34.548 

20.725 

376.000 

34.865 

21.083 

380.000 

35.188 

21.429 

384.000 

35.536 

21.786 

388.000 

35.916 

22.138 

392.000 

36.254 

22.503 

396.000 

36.571 

22.860 

400.000 

36.881 

23.182 

404.000 

37.206 

23.534 

408.000 

37.547 

23.905 

412.000 

37.886 

24.270 

416.000 

38.223 

24.606 

420.000 

38.558 

24.937 

424.000 

38.892 

25.289 

428.000 

39.225 

25.652 

432.000 

39.557 

25.989 

436.000 

39.889 

26.337 

440.000 

40.219 

26.693 

444.000 

40.549 

27.026 

448.000 

40.899 

27.361 

452.000 

40,706 

27.715 

456.000 

40.098 

28.075 

460.000 

42.645 

28.407 

464.000 

42.382 

28.748 

468.000 

42.644 

29.099 

472.000 

42.908 

29.438 

476.000 

43.180 

29.575 
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43.566 

29.897 

484.000 

43.809 

30.714 

488.000 

44.120 

30.943 

492.000 

44.355 
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RANGE 


PLOW 

PHIGH 

COEF 

EXP 

CORR  COEF 

MPA 

MPA 

CM/S-MPA^EXP 

— 

25. 

50. 

.164721E+00 

.901040 

.999392 

50. 

124. 

.160656E+00 

.907147 

.999992 

124. 

372. 

.159972E+00 

.908079 

.999962 

10. 

490. 

.160746E+00 

.907204 

.999988 

125 


iNTEKnONALLY  LEFT  BLANK. 


.. 

• 

No.  of 

No.  of 

Conies 

Orcani7.ation 

Conies 

Organization 

2 

Administrator 

1 

Commander 

Defense  Technical  Info  Center 

U.S.  Army  Missile  Command 

ATTN:  DTIC-DDA 

ATTN:  AMSMI-RD-CS-R  (DOC) 

■  : 

Cameron  Station 

Redstone  .Arsenal,  AL  35898-5010 

Alexandria,  VA  22304-6145 

1 

Commander 

1 

Commander 

U.S.  Army  Tank-Automotive  Command 

. 

U.S.  Army  Materiel  Command 

ATTN:  ASQNC-TAC-DIT  (Technical 

ATTN:  AMCAM 

Information  Center) 

5001  Eisenhower  Ave. 

Warren,  MI  48397-5000 

Alexandria.  VA  22333-0001 

i 

1 

Director 

1 

Director 

U.S.  Army  TRADOC  Analysis  Command 

U.S.  Army  Research  Laboratory 

ATTN:  ATRC-WSR 

ATTN:  AMSRL-D 

White  Sands  Missile  Range,  NM  88002-5502 

2800  Powder  Mill  Rd. 

1 

Adelphi,  MD  20783-1145 

1 

Commandant 

U.S.  Army  Field  Artillery  School 

1 

Dircctoi 

ATTN:  ATSF-CSI 

.  5  ■  • 

U.S.  Army  Research  Laboratory 

Ft.  SUl,OK  73503-5000 

ATTN:  AMSRL-OP-CI-AD, 

Tech  Publishing 

(CliM.  only)  1 

Commandant 

2800  Powder  fvIUl  Rd. 

U.S.  Anny  Infantry  School 

Adelphi,  MD  20783-1145 

ATTN:  ATSH-CD  (Security  Mgr.) 

Fort  Denning,  GA  31905-5660 

2 

Commander 

U.S.  A-my  Armament  Research, 

(UncItM.  only)  ] 

Commandant 

Development,  and  Engineering  Center 

U.S.  Army  Infantry  School 

ATTN:  SMCAR-IMI-I 

ATTN:  ATSH-CD-CSO-OR 

V  1 

Picatinny  Arsenal,  NJ  07806-5000 

Fort  Denning,  GA  31905-5660 

( 

1 

2 

Commander 

1 

WL/MNOI 

U.S.  Army  Armament  Research, 

Eg’in  AFD,  FL  32542-5000 

•41 

Development,  and  Engineering  Center 

A'fTN:  SMCAR-TDC 

Aberdeen  Proving  Ground 

Picatinny  Arsenal,  NJ  07806-5000 

2 

Dir,  USAMSAA 

\ 

Director 

ATTN;  AMXSY-D 

Benet  Weapons  Laboratory 

AMXSY-MP.  H.  Cohen 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 

1 

Cdr,  USATTCOM 

ATTN:  SMCAR-CCB-TL 

ATTN:  AMST1E  TC 

..  i 

Watcrvlict,  NY  12189-4050 

1 

Dir,  ERDEC 

;■  ■ 

(Unclatt.  only)  | 

Commander 

ATTKJ* 

U.S.  Army  Rock  Island  Asenal 

ATTN:  SMCRI-IMC-RT/Tcchnical  Library 

1 

Cdr,  CBDA 

Rock  Island,  IL  61299-5000 

ATTN:  AMSCB-CI 

1 

Director 

1 

Dir,  USARL 

U.S.  Army  Aviation  Research 

ATTN:  AMSRL-SL-I 

1 

and  Technology  Activity 

1 

ATTN:  SAVRT-R  (Library) 

10 

Dir,  USARL 

1 

M,/S  219-3 

ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 

’ 

Ames  Reseuch  Center 

" 

Moffett  Field,  CA  94035-1000 

1. 
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1  OSD/SDIOAST 

ATTN:  Dr.  Len  Caveny 
Pentagon 

Washington.  DC  20301-7100 
1  Commander 

U.S.  Army  Aimamcnt  Research,  Development, 
and  Engineering  Center 
ATIN:  SMCAR-QAH-T,  Mi.  John  Domcn 
Bldg.  62  North 

Picatinny  Arsenal,  NJ  07805-5000 
1  Commander 

U.S.  Army  Armament  Research,  Development, 
and  Engineering  Center 
ATTN:  SMCAR-AEE-BR, 

Mr.  Ken  Klingaman 

Bldg.  1501 

Picatinny  Arsenal,  NJ  07806-5000 
1  Commander 

U.S.  Army  Armament  Research,  Development, 
and  Engineering  Center 
ATTN:  SMCAR-ASI,  Dr.  Jerry  Rubin 
Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

U.S.  Army  BMD 
Advanced  Technology  Center 
P.  O.  Box  1500 
Huntsville,  AL  35807 

1  Chairman 

DOD  Explosives  Safety  Board 
Room  856<; 

Hoffman  Bldg.  1 
2461  Eisenhower  Ave. 

Alexandria.  V A  22331-0600 

1  Department  of  the  Army 

Office  of  the  Product  Manager 
155mm  Howitzer,  M109A6,  Paladin 
AITN:  SFAE-AR-HIP-IP,  Mr.  R.  Dc  Kleinc 
Picatinny  Aiscnai,  NJ  07806-5000 


1  Commander 

Production  Base  Modeniization  Agency 
U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  AMSMC-PBM-E.  L.  Laibson 
Picatinny  Arsenal,  NJ  07806-5000 

3  PEO-Armaments 
Project  Manager 
Tank  Main  Armament  Systems 
ATTN:  AMCPM-TMA,  K.  Russell 
AMCPM-TMA-105 
AMCPM-TMA- 120 
Picatinny  Arsenal,  NJ  07806-5000 

6  Director 

Benet  Laboratories 

U.S.  Army  Watcrviict  Arsenal 

ATTN:  SARWV-RD, 

L.  Johnson 
G.  Carafano 
R.  Thierry 
R.  Hasoenbein 
P.  Vous 
P.  Alto 

Watcrviict,  NY  12189 

3  Commander 

U.S.  Anny  AMCCOM 
ATTN:  AMSMC-IRC,  G.  Cowan 
SMCAR-ESM(R), 

W.  Fortune 
R.  Zastrow 

Rock  Island,  IL  61299-7300 

1  Cosnmander,  USACECOM 
R&D  Technical  Library 
ATTN:  ASQNC-ELC-IS-L-R,  Myer  Center 
Fort  Monmouth,  NJ  07703-5301 

1  Commandant 

U.S.  Army  Aviation  School 
ATTN:  Aviation  Agency 
Fort  Rucker,  AL  36360 

1  Director 

HQ,  TRAC  RPD 

ATTN’:  ATCD-MA,  MAJ  Williams 
Fort  Monroe,  VA  23651-5143 
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1  Headquarters 

U.S.  Army  Materiel  Command 
ATTN:  AMCICP-AD,  Michael  F.  Fisette 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

4  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCD,  D.  Spring 
SMCAR-CCS 

SMCAR-CCH-T,  L.  Roscndorf 
SMCAR-CCH-V,  E.  Fennell 
Picatinny  Arsenal,  NJ  07806-5000 

8  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-AE,  J.  Picard 
SMCAR-AEE-B, 

A.  Bcardell 
D.  Downs 

S.  Einstein 

A.  Bracuti 
D.  Chiu 

SMCAR-AEE,  J.  Lannon 
SMCAR-AES,  S.  Kaplowitz 
Picatinny  Arsenal,  NJ  07806-5000 

7  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-FSA-T,  M.  Salsbury 
SMCAR-FSE, 

T.  Gora 

B.  Knutelsky 
K.  C.  Pan 
W.  Davis 

C.  Durham 
A.  Graf 

Picatinny  Arsenal,  NJ  07806-5000 
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Copies  Organization 

6  Commairder 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN:  SMCAR-EG, 

G.  Ferdinand 

H.  Nabcr-Libby 
R.  Lundberg 

J,  Niles 
R.  Moteira 
W.  Morclli  • 

Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

U.S.  Army  Research  Office 
ATTN:  Technical  Library 

D.  Mann 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-221 1 

1  Commander 

U.S.  Arm>  Bclvoir  R&D  Center 
ATTN:  STRBE-WC, 

Technical  Library  (Vault) 

Bldg.  315 

Fort  Bclvoir,  VA  22060-5606 

1  Commander 

U.S.  Anny  TRAC  •  Foit  Lee 
Defense  Logistics  Studies 
Fort  Lee,  VA  23801-6140 

1  Pi'csidcnt 

U.S.  Army  Artillery  Board 
Fort  Sill,  OK  73503 

1  Commandant 

U.S.  Army  Command  and  General  Staff  College 
Fort  Leavenworth,  KS  66027-5200 

1  Commandant 

U.S.  Army  Special  Warfare  School 
ATTN:  Rev  and  Tng  Lit  Div 
Fort  Bragg,  NC  28307 

1  Commander 

Radford  Army  Ammunition  Plant 
ATTN:  SMCRA-QA,  HI  Librarj' 

Radford.  VA  24141 
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1  Commandant 

U.S.  Army  Field  Artillery  School 
ATTN:  STSF-TSM-CN 
Fort  Sill,  OK  73503-5600 

4  Deputy  Commander 

Strategic  Defense  Command 
ATTN:  SFAE-SD-HVL, 

S.  Smith 
LTC  Kce 
D.  Lianos 

T.  Aden 
P.O.  Box  1500 
Huntsville.  AL  35887-8801 

3  Commander 

U.S.  Army  Foreign  Science  and  Technology 
Center 

ATTN:  AMXST-MC-3. 

S.  LcBeau 
C.  Beitcr 
220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

1  Commandant 

U.S.  Army  Field  Artillery  Center  and  School 
ATTN:  ATSF-CO-MW,  B.  Willis 
FortSUl.OK  73503 

1  Naval  Sea  System  Command 
Dept  of  the  Navy 
AT^:  CSEA,  CDR  Dampier 
OdKRr- 

Washington,  DC  20362-5101 

1  Office  of  Naval  Research 
AITN:  Code  473,  R.  S.  Miller 
800  N.  Quincy  St 
Ailingion,  VA  22217 

2  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA62R 
SEA  64 

Washington,  DC  20362-5101 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-954,  Technical  Library 
Washington,  DC  20360 


1  Naval  Research  Laboratory 
Technical  Library 
Washington,  DC  20375 

2  Commander 

Naval  Surface  Warfare  Center 
ATTN:  J.  P.  Consaga 
C.  Gotzmer 

Silver  Spring,  MD  20902-5000 

2  Commander 

Naval  Surface  Warfare  Center 
ATTN:  Code  R-13.  K.  Kim 

Code  R-13,  R.  Bcmcckcr 
Silver  Spring,  MD  20902-5000 

3  Commander 

Naval  Surface  Warfare  Center 
Indian  Head  Division 
ATTN:  610,  C.  Smith 
6110J,  K.  Rice 
61  IOC,  S.  Peters 
Indian  Head,  MD  20640-5035 

6  Commander 

Naval  Surface  Warfare  Center 
Dahlgicn  Division 
ATTN:  CodcG33, 

T.  Doran 
J.  Copley 

Code  G30,  Guns  and  Munitions  Division 
Code  G301.  D.  Wilson 
Code  G32,  Guns  Systems  Division 
Code  E23,  Technical  Library 
Dahlgrcn,  VA  22448-5000 

1  Commander 

Naval  Surface  Warfare  Center 
Crane  Division 

ATTN:  Code  4052,  Steve  Backer 
Bldg.  108 

Crane,  IN  47522-.5000 

2  Commander 

Naval  Underwater  Systems  Center 
Energy  Convention  Dept 
ATTN:  Code  5B331.  R.  S.  Lazar 
Technical  Library 
Newport,  RI  02840 
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1  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
ATTN;  Codc270Pl.  Mr.  Ed  Chan 
10 1  Straus  Avc. 

Indian  Head,  MD  20640 

1  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
ATTN:  Code  3120,  Mr.  Robert  Rast 
101  Straus  Avc. 

Indian  Head,  MD  20640 

1  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
ATTN:  Code  210P1,  Mr.  Ron  Simmons 
101  Straus  Avc. 

Indian  Head,  MD  20640 

2  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
A7TN:  Code  6210, 

Sharon  Boyles 
Norberto  Almcyda 
101  Straus  Avc. 

Indian  Head,  MD  20640 

1  Commatidcr 

Naval  Air  Warfare  Center 
ATIX:  Code  3891,  Mr.  Chan  Price 
China  Lake,  CA  93555 

1  Commander 

Naval  j\ir  Warfare  Center 

ATTN:  Code  3891,  Ms.  Alice  Atwood 

China  Lake,  CA  93555 

1  Commander 

Naval  Weapons  Center 
ATTN:  Code  388,  C.  F.  Price 
Info  Science  Division 
China  Lake,  CA  93555-6001 

1  OLAC  PL/rSTL 
ATTN:  D.  Shiplctt 
Edwards  AI^,  CA  93523  -5000 


No.  of 

Conies  Organiration 

10  Central  Intelligence  Agency 

Office  of  Central  Reference  Dissemination 
Branch 

Rm.  GE-47  HQS 
Washington,  DC  20502 

1  Cenual  Intelligence  Agency 

ATTN:  Joseph  E.  Backofen 
HQ  Rm.  5F22 
Washington,  DC  20505 

5  Director 

Sandia  National  Laboratories 
ATTN:  T.  Hitchcock 

R.  Woodfin 
D.  Benson 

S.  Kempka 

R.  Beasley 

Advanced  Projects  Division  14 
Organization  9123 
Albuquerque,  NM  87185 

2  Director 

Los  Alamos  National  Laboratory 
ATTN:  B.  Kaswhia 
H.  Davis 

Los  Alamos,  NM  87545 

1  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  M.  S.  L-355,  A.  Buckingham 
P.O.  Box  808 
Livermore,  CA  94530 

2  Dirccujr 

Sandia  National  Laboratories 
Combustion  Research  Facility 
ATTN:  R.  Armstrong 

S.  Vosen 
Division  8357 

Livermore,  CA  94551-0469 

1  University  of  Illinois 

Dept  of  Mcch./IndusU  Engr. 

ATTN:  Professor  Herman  Kricr,  144  MEB 
1206  N.  Green  St. 

Urbana,  IL  61801 
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1  The  Johns  Hopkins  University/CPIA 
ATTN:  T.  Christian 

10630  Little  Patuxent  Parkway,  Suite  202 
Columbia,  MD  21044-3200 

2  Pennsylvania  State  University 
DepL  of  Mechanical  Engr. 

ATTN:  Jeff  Brown 

312  Mechanical  Engineering  Bldg, 
University  Park,  PA  16802 

1  Noith  Carolina  State  University 
ATTN:  John  G.  Gilligan 

Box  7909 

1110  Burlington  Engineering  Labs 
Raleigh,  NC  27695-7909 

2  Institute  for  Advanced  Studies 
ATTN:  Dr.  H.  Fair 

Dr.  T.  Kiehne 
4030-2  West  Baker  Lane 
AusUn,  TX  78759-5329 

1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Technical  Library 
333  Ravenswood  Avc. 

Menlo  Park,  CA  94025 

1  SPARTA 

ATTN:  Dr.  Michael  Holland 
9455  Towne  Center  Dr. 

San  Diego,  CA  92121-1964 

5  FMC  Corporation 

ATTN:  Mr.  G.  Johnson 
Mr.  M.  Scale 
Dr.  A.  GiovanctU 
Mr.  J.  Dyvik 
Dr.  D.  Cook 
4800  East  River  Rd. 

Minneapolis,  MN  55421-1498 


2  Hercules  Inc. 

Radford  Army  Ammunition  Plant 
Manager  Manufactuiing  Engineering  Dept. 
ATTN:  D.  A.  WorreU 
Edward  Sanford 
P.O.  Box  1 
Radford,  VA  24141 

1  Hercules  Inc. 

ATTN:  Dr.  Richard  Cartwright 
100  Howard  Blvd. 

Kenvil.NJ  07847 

3  GT  Devices 

ATTN:  Dr.  S.  Goldstein 
Dr.  R.  J.  Grcig 
Dr.  N.  Winsor 

5705A  General  Washington  Dr. 
Alexandria,  VA  22312 

3  General  Dynamics  Land  Systems 
ATTN:  Dr.  B.  VanDcusen 

Mr.  F.  Lunsford 
Dr.  M.  Weidner 
P.O.  Box  2074 
Warren,  Ml  48090-2074 

2  Alliant  Tcchsystcms,  Inc. 

ATTN:  R  .  E.  Tompkins 

J.  Kennedy 
7225  Northland  Dr. 

Brooklyn  Park,  MN  55428 

4  Olui  Ordnance 

ATTN:  V.  McDonald,  Library 
Hugh  McElroy 
Mr.  Tiiomas  Bourgeois 
Mr .  Dennis  Wortliington 
P.O.  Box  222 
St.  Marks,  FL  32355 

1  Paul  Gough  Associates,  Inc. 

ATTN:  P.  S.  Gough 
1048  South  St. 

Portsmouth,  NH  03801-5423 
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1  Physics  Intcmatiotial  Library 
ATTN:  H.  Wayne  Wampler 
P.O.  Box  5010 

San  Leandro,  CA  94577-0599 


1  Veritay  Technology,  Inc. 
ATTN;  Ivfr.  E.  Fisher 
4845  Millersport  Highway 
East  Amherst,  NY  14051-0305 


2  Rockwell  International 
Rocketdyne  Division 
ATTN:  BA08, 

J.  E.  Flanagan 
J.  Gray 
6633  Canoga  Ave. 

Canoga  Park,  CA  91304 

2  Princeton  Combustion  Research 
Laboratories,  Inc. 

ATTN:  M.  SummeiTield 
N.  Messina 

Princeton  Corporate  Plaza 
1 1  Deerpark  Dr. 

Bldg.  IV,  SuiU*  1 19 
Monmouth  Junction,  NJ  08852 


1  Veritay  Technology,  Inc. 

4845  Ivlillcrsport  Highway 
P.O.  Box  305 

East  Amherst,  NY  14051-0305 

1  Baitcllc 

ATTN:  TACTEC  Library,  J.  N.  Huggins 
505  IGng  Ave. 

Columbus,  OH  43201-2693 

2  California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
ATTN:  L.  D.  Strand,  MS  125-224 

D.  Elliot 

4800  Oak  Grove  Dr. 

Pasade  la,  CA  91109 


2  Science  Applications,  Inc. 
ATTN:  J.  Battch 

L.  Thornhill 
1519  Johnson  Ferry  Rd. 
Suite  300 

Marietta,  GA  30062-6438 


1  General  Electric  Co. 
Dcfcn.se  Systems  Division 
ATTN:  Dr.  J.  Mandzy 
Mail  Drop  43-220 
100  Plastics  Ave. 
Pittsfield,  MA  01201 


1  Eli  Freedman  &  Associates 
ATTN:  E.  Freedman 
2411  Diana  Rd. 

Baltimore,  MD  21209 

1  Rocketdyne 

ATTN:  Mr.  Otto  Hcincy 
Mail  Stop  BA26 
6633  Canoga  Ave. 

Conoga  Park,  CA  91304 

1  Thiokol  (Longhorn  Division) 

AITN:  Dr.  David  Dillchay 
Mail  Stop  703-1 1 
P.O.  Box  1 149 
Marshall,  TX  75671 

1  Thiokol  (Elkton  Division) 

ATTN:  Dr.  Rodney  Wilier 
55  Thiokol  Rd. 

Elkton,  MD  21922 


2  SAIC 

ATTN:  Mr.  N.  Sinha 
Dr.  S.  Dash 
501  Office  Center  Dr. 

Fort  Washington,  PA  19034-3211 

1  State  University  of  New  York 

Department  of  Electrical  Engineering 
ATTN:  Dr.  W.  J.  Sargeant 
Bonner  Hall  -  Rm.  312 
Buffalo,  NY  14260 

1  Science  Application  Intematicnal  Corporation 

ATTN:  Dr.  George  Chryssamdlis 
8400  Normandelc  Blvd. 

Suite  939 

Minneapolis,  MN  55437 
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1  Hercules,  Inc. 

ATTN:  Dr.  Richard  Cartwright 
100  Howard  Blvd. 

Kcnvil.NJ  07847 

1  Sandia  National  Laboratories 

ATTN:  Mr.  Mark  GrubcUch,  DIV  2515 
P.O.  Box  5800 
Albuquerque,  NM  87185 

Abt‘-rdeen  Proving  Ground 

4  Cdr,  USACSTA 
ATTN:  S.  Walton 
G.  Rice 
■  D.  Lacey 
C.  Hcrud 


No.  of 
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2  RARDE 

GS2  Division 
Bldg.  R31 

ATTN:  Dr.  C.  Woodley 
Dr.  G.  Cook 
Fort  Halstead 

Sevenoaks,  Kent  TN14  7BP 
ENGLAND 

1  Materials  Research  Laboratory 

Salisbury  Branch 

ATTN:  Anna  Wildcgger  Gaissinaicr 
Explosives  Ordnance  Division 
Salisbury 

SOUTH  AUSTOALIA  5108 

1  Laboratorio  Quimico  Central  de  Armaincnto 

AITN:  Captain  Juan  F.  Hernandez  Tamayo 
Apartado  1105 
28080  Madrid 
SPAIN 

1  R&D  Dept. 

ATTN:  Dr.  Pierre  Archambault 
5  Montec  dcs  Arsenaux 
Le  Gardeur,  Quebec 
CANADA  J5Z2P4 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


Ttiis  Laborator)'  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-36  (Part  I) _ Date  of  Report  January  1993 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need  ■  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 

4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 

ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS _ _ _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
abovt  and  the  Old  or  Incorrect  address  below. 


Organization 

OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


